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ABSTRACT 

We have spectroscopically identified 60 Galactic Wolf-Rayet stars (WRs), including 38 nitrogen types 
(WN) and 22 carbon types (WC). Using photometry from the Spiteer /GLIMPSE and Two-Micron 
All-Sky Survey (2MASS) databases, the new WRs were selected via a method we have established 
that exploits their unique infrared colors, which is mainly the result of excess radiation generated by 
free-free scattering within their dense ionized winds. The selection criterion has been refined since 
the last report, resulting in a WR detection rate of w20% in spectroscopic follow-up of candidates 
that comprise a broad color space defined by the color distribution of all known WRs having B > 14 
mag. However, there are smaller regions within this color space which yield WRs at a rate of >50% 
in spectroscopic follow-up. Candidates which are not WRs are mainly Be stars, which is possibly 
attributable to the physical similarities between the free-free emission parameters of Be disks and 
WR winds. As an additional selection experiment, the list of WR candidates was cross-correlated 
with archival X-ray point-source catalogs, which increases the WR detection rate of the broad color 
space to ten new WR X-ray sources have been found, in addition to a previously unrecognized 

X-ray counterpart to a known WR. The extinction values, distances, and galactocentric radii of all 
new WRs are calculated using the method of spectroscopic parallax. Although the majority of the 
new WRs have no obvious association with stellar clusters, two WC8 stars reside in a previously 
unknown massive-star cluster, in which five OB supergiants were also identified. The new system lies 
at an estimated distance of « 6.1 kpc, near the intersection of the Scutum-Ccntaurus Arm with the 
Galaxy's bar. In addition, two WC and four WN stars, all but one of which are X-ray sources, were 
identified in association with the stellar clusters Danks 1 and 2. A WN9 star has also been associated 
with the cluster [DBS2003] 179. This work brings the total number of known Galactic WRs to 476, or 
«7%-8% of the total empirically estimated population. An examination of their Galactic distribution 
reveals an approximate tracing of spiral arms and an enhanced WR surface density toward several 
massive-star formation sites. 
Subject headings: 



1. INTRODUCTION 

Wolf-Rayet (WR) stars are the evolved descendants 
of the most massive stars (M > 25 M©; see Crowther 
2007 and references therein). They are highly luminous 
{L ~ lO^-lO^i©) and drive ferocious winds, which, when 
integrated over the entire WR phase, can deposit an 
amount of mechanical energy comparable to that released 
by a supernova blast. As the strong winds of WRs eject 
hydrogen and helium from their atmospheres, deeper 
and hotter layers containing the byproducts of core nu- 
clear burning arc exposed. The high temperatures and 
fast winds result in the development of strong, velocity- 
broadened emission lines of helium, carbon, nitrogen, 
and oxygen. Based on the observed relative strengths of 
these lines, WRs are classified into spectral subtype cat- 
egories of nitrogen type (WN), carbon type (WC), and 
oxygen type (WO). These subtypes are believed to fol- 
low the general evolutionary sequence: WN^WC— t-WO 
(e.g., Langer et al. 1994; Norci et al. 1998). As WRs are 
hydrogen deficient upon death, they are suspected to be 
the progenitors of Type lb and Ic core-collapse super- 
novae and long-duration gamma-ray bursts. With ages 
of <10 Myr, WRs are excellent markers of recent massive 
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star formation. Their population demographics provide 
insight on the formation mode and distribution of the 
most massive stars in the Galaxy, while the observed fre- 
quencies of WR subtypes as a function of environmental 
metallicity provide a useful testing ground for evolution- 
ary models (Meynet & Maeder 2005). 

Prior to this work, there were upward of 400 cata- 
loged WRs. Those that lie relatively close to the Sun 
have bright optical counterparts (van der Hucht 2001) 
and this local population is estimated to be (;omplete to 
B < lA mag (Shara et al. 1999). Based on the local 
WR space density, a simple extrapolation implies that 
«6500 WRs populate the Galaxy (van der Hucht 2001; 
Shara et al. 1999, 2009). Since the hfetime of WR precur- 
sors is <10 Myr, they are unable to drift far from their 
birthplaces in molecular clouds before the end of their 
lives. Hence, from the vantage point of Earth, the vast 
majority of Galactic WRs lie in regions of high optical 
obscuration. Fortunately, infrared observations have the 
ability to unveil WRs in such environments. Large popu- 
lations of obscured WRs have been revealed in regions of 
the Galaxy for which astronomical interest is relatively 
biased, such as the Galactic center (Figer et al. 1999, 
2002; Mauerhan et al. 2007, 2009, 2010b) and extraordi- 
nary starburst clusters such as Westerlund 1 (Crowther 
et al. 2006). More recently, infrared-based searches for 
WRs have been extended to the greater Galactic plane. 
Shara et al. (2009) have identified 41 obscured WRs 
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via infrared narrow-band imaging and follow-up spec- 
troscopy. Meanwhile, complementary efforts by Hadficld 
et al. (2007) and Mauerhan at al. (2009, 2010a) have 
utilized the 2MASS and Spitzer/ GLIMPSE point-source 
catalogs to develop a very efficient color-based criterion 
for selecting WR candidates from the field. This work is 
a continuation of that effort. We report the detection of 
60 new WRs, which brings the current total number of 
WRs selected by our method to ssQO. The color-based 
criterion has been refined since our last major discovery 
report in Mauerhan et al. (2009), leading to an enhanced 
detection rate for WRs and a better suppression of non- 
WR "contaminating" sources. 

2. WR SELECTION CRITERION 
2.1. Infrared Color and Brightness 

WRs are significant sources of infrared excess. Figure 
1 demonstrates that they are clearly separable from the 
vast majority of field stars in {H — Kg) versus (J — H), 
{Kg - [8.0]) versus (J - Kg), and ([3.6] - [8.0]) versus 
([3.6] — [4.5]) color-color diagrams. The color excesses 
are mainly the result of free-free emission from e~ scat- 
tering off of helium ions in the outer regions of their 
dense winds (Wright & Barlow 1975; Cohen et al. 1975). 
Broad emission lines and thermal emission from hot dust 
for a subset of WC stars can provide additional color ex- 
cess (e.g., see Williams et al. 1987). To illustrate the 
extent by which line emission contributes to the infrared 
excess of WRs, we removed the emission lines from the 
flux-calibrated spectra of selected objects in our sam- 
ple and used the STSDAS SYNPHOT package in IRAF^ 
to generate synthetic JHKg photometry from the pure- 
continuum photometry. Figure 2 illustrates the resulting 
shift in {H — Kg) and (J — H) color that several WRs ex- 
perience after their line emission was removed. Selected 
WN stars from our sample do not undergo a significant 
color shift as a result of emission lines, which probably 
explains why WN stars form a relatively tight trend in 
the {H — Kg) versus (J — H) color space of Figure 1. 
However, WC stars undergo a significant shift m. H — Kg 
color if their emission lines are removed. The effect is 
mainly driven by very strong emission from C IV and C 
III in the /C-band. Relative to the continuum, WR line 
strengths decrease with increasing wavelength, so, line 
emission has a much smaller eflFect on broad-band colors 
of WRs in the mid-infrared. 

The infrared excesses of WRs shown forms the basis 
of a selection criterion within the area of sky covered 
by the 2MASS and 5piteer/GLIMPSE surveys of the in- 
ner Galaxy [l = 295°-0°, 0°-65°; h = ±1°; Benjamin 
et al. 2003). The color space used in this work has 
been expanded to encompass more distant and highly 
obscured WR candidates that lie along the reddening 
vector. Furthermore, we determined that the inclu- 
sion of color-magnitude restrictions provides an addi- 
tional "sieve" for separating WRs from non-WR "con- 
taminants". The color spaces bounded by solid lines in 
Figure 1 were defined by the color distributions of all 
known WRs having _£? > 14 mag (the Galactic WR sam- 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which arc operated by the Association of Universities 
for Research in Astronomy, Inc. , under cooperative agreement with 
the National Science Foundation. 



pie is already complete to i? < 14 mag; Shara et al. 1999) 
We will hereafter refer to the spaces defined by the black 
lines in Figure 1, collectively, as the broad color space. 
Known WRs of all subtypes occupy this space, including 
dusty WC stars. There are ^6600 unidentified WR can- 
didates in the space having Kg < 13.1 mag; these have 
become targets for spectroscopic characterization. Fo- 
cusing on the brightest sources with K < 11.5 mag, we 
have experienced a Ri20% detection rate for WRs in the 
broad color space. In an attempt to increase the WR de- 
tection rate, we experimented with selecting candidates 
from smaller color subspaces in the (Kg — [8.0]) versus 
(J - Kg), and ([3.6] - [8.0]) versus ([3.6] - [4.5]) spaces 
where non-WR "contaminants" are relatively sparse; we 
will refer to these color subspaces, collectively, as the 
"VFi? sweet spot" . As expected, the WR detection rate 
increased significantly to upward of s=;50%, and the rate 
of WC type detection also increased notably. 

It is not entirely clear what infrared spectral and ex- 
tinction properties drive WR colors into the sweet spot or 
what subtype selection biases result from selecting candi- 
dates from this color region. A clue might lie in the fact 
that the detection rate for WC stars in the sweet spot 
increased relative to the detection rate using the broad 
color criterion. The strong line emission from C III and 
C IV in the JsTs-band could explain the sweet spot in the 
(iiTs — [8.0]) versus (J — Kg) color space (moving strong- 
lined WC stars blueward in iirs-[8.0] axis and redward 
in J — Kg axis), but emission- line effects are probably 
less significant for the [3.6] — [8.0] versus [3.6] — [4.5], 
because the line strengths of WRs relative to the contin- 
uua decreases significantly in the mid-infrared. It is also 
plausible that the WR sweet spot is produced, in part, by 
the relative strength of free-free emission for particular 
subtypes, which could be the cause of the enhanced de- 
tection rate of WC types in this subspace. Whatever the 
case, it should be noted that representatives from most 
WN and WC subtypes are found within the WR sweet 
spot, which suggests that there is no exclusive subtype 
bias in this particular color space with respect to the 
broader color space. However, the color distribution of 
dusty WC stars in Figure 1 suggests that there will be 
some bias against these particular types in the WR sweet 
spot. 

Additional work is needed to establish an improved 
understanding of the WR phenomenology in the infrared 
and the selection biases and limitations of our survey. 
The main goal of this work was to maximize the detection 
rate for Galactic WRs, in order to best exploit the limited 
observing time made available for our program. We have, 
thus far, given precedence to the brightest WRs candi- 
dates in our sample (K < 11.5 mag), in order to expedite 
the discovery of new WRs. As a result, we have limited 
the Galactic volume probed by our survey. For instance, 
by focusing the spectroscopic sample to objects having 
Kg < 11.5 mag, we are mainly sensitive to WRs lying 
at distances of <8 kpc, assuming an average Ak^ = 2 
and an average absolute magnitude of M^^ = —5 for 
WRs (see Crowther et al. 2006). According to Shara et 
al. (2009), «2500 Galactic WRs (^40%) are expected to 
have K > 13 mag. Most of these fainter WRs, which we 
have not yet sampled in our survey, will lie on the far 
half of the Galaxy. Thus our survey is strongly biased 
towards WR candidates on the near side of the Galaxy. 
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Most of the non-WR stars selected in our survey ex- 
hibit the spectral characteristic of Be stars. These ob- 
jects exhibit H I Paschen and Brackett series emission in 
their spectra, shown in Figure A.l, which is consistent 
with the H and K spectral atlases of Be III-V stars from 
Clark & Steele (2000) and Steele & Clark (2001). We sus- 
pect the reason that WRs and Be stars have such similar 
infrared colors is a result of the similarities between the 
sources of their infrared excesses, which, in the case of Be 
stars, is free-free emission generated by their ionized cir- 
cumstellar disks. For both WR winds and Be disks, the 
frcc-frcc flux is a function of electron density (iVg), elec- 
tron temperature (T^), the physical size of the emitting 
source {Rs), and the wavelength at which the emission 
becomes optically thick (Aj). All of these parameters 
are practically identical for WRs (see Cohen et al. 1975; 
Mauerhan et al. 2010b) and Be stars (see Gchrz ct al. 
1974), having typical values on the order of A^e ^ 10^^ 
cm-3, Te - lO'^ K, Rtt = 10^2-1013 cm, and A* « 5 /im. 
These physical similarities provide a simple and natural 
explanation for the infrared color spaces that WRs and 
Be stars share. In addition to Be stars, our method has 
also selected several candidate evolved B[e] and O su- 
pergiants, and we suspect that these objects were also 
selected because of similarities in the free-free parame- 
ters. The discovery of these stars will be presented in a 
forthcoming paper. 

The refinement of the selection method in our most 
recent survey resulted in a significant improvement over 
our previous WR return rate of <10% reported in Mauer- 
han et al. (2009). It appears that by including the recent 
color-magnitude restriction, the frequency of cool, late- 
type "contaminants" has dropped significantly (to < 5%) 
since the surveys of Hadfield et al. (2007) and Mauer- 
han et al. (2009). The broad color space, thus, returns 
early-type emission-line stars almost entirely, at a rate of 
«95%. 

2.2. X-ray Emission: An Additional Selection Tool 

WRs are commonly X-ray sources (see Skinner et al. 
2010 and references therein). The high-energy emission 

is generated via shocks that form from line-driven in- 
stabilities within their winds (Lucy & White 1980), or 
from the collision between opposing winds in a massive 
binary (Luo et al. 1990). Thus, as a complementary 
WR selection technique, we cross-correlated our list of 
infrared color-selected WR candidates from the broad 
color space of Figure 1 with archival X-ray catalogs that 
were compiled from multiple Galactic plane observations 
made with the Chandra X-ray Observatory (Evans et al. 
2010) and XMM-Nemton (Hands et al. 2004). We recov- 
ered 43 WR-candidate/X-ray matches, including several 
known WRs that were not prcvioiisly recognized as X-ray 
sources. We obtained spectra for 23 of the candidates and 
recovered 10 WRs, for a return rate of ~40%, a factor-of- 
two increase in the return rate using the broad infrared 
color space alone. The non-WRs in this particular sam- 
ple are, still, mainly Be stars, but we also identified sev- 
eral candidate YSOs, and late-type stars with infrared 
excess; these objects will be the subject of a forthcoming 
paper. The basic X-ray properties of these new WRs are 
presented in Section 6. 

Although X-ray emission proves to be an effective tool 
for selecting WRs from the list of infrared-color candi- 



dates and improving the WR detection rate, there are 
certainly biases towards certain WR subtypes. For in- 
stance, although X-ray emission appears to be ubiqui- 
tous property of WN stars (Skinner et al. 2010), most 
WC stars are undetectable as X-ray sources (Oskinova 
et al. 2003), unless they are members of colliding- wind 
binary systems. Furthermore, the X-ray emission from 
single WRs is typically dominated by a soft X-ray pho- 
tons having energies below 1 keV, which are heavily ab- 
sorbed by the interstellar medium. Thus, X-ray selection 
will be biased toward WRs in regions of lower extinction, 
or those that emit harder X-rays (e.g., colliding- wind bi- 
naries) which are less attenuated by dust and gas along 
the intervening ISM. 

3. INFRARED SPECTROSCOPY 

Near-infrared spectra of WR candidates was obtained 
with the Palomar Hale 5 m telescope and the TripleSpec 
(TSPEC) near-infrared spectrograph during 2009 July 
13-15 and 2010 June 21 (UT), the NASA Infrared Tele- 
scope Facility (IRTF) and SpeX instrument during 2009 
August 4 6 (UT). and the Southern Observatory for As- 
trophysical Research (SOAR) and OSIRIS Spectrograph 
during 2010 May 23 27 (UT). TSPEC provides simul- 
taneous coverage of the J, H, and K bands (A = 1.0- 
2.4 /im) and produces a moderate-resolution spectrum 
(i? w 2500-2700) through a 1" x 30" slit. SpcX and 
OSIRIS provide wavelength coverage similar to TSPEC, 
although at respective resolutions of i? w 1200 through a 
0'.'5 slit and R « 800 through a 1" slit. All spectra were 
obtained utilizing an ABBA telescope nodding sequence 
for the subtraction of sky emission and to suppress bad 
pixels. For the TSPEC data, wavelength calibration was 
performed using the OH emission lines detected in the 
background sky spectra. For the SpeX and OSIRIS data, 
spectra from arc lamps internal to the instruments were 
used for wavelength calibration. TSPEC and OSIRIS im- 
ages were reduced using standard IRAF routines, and the 
spectra were extracted using the IRAF package APALL. 
The SpeX data were reduced using the IDL-based pack- 
age Spextool (Cushing et al. 2004). For all observations 
AOV standard stars were observed to obtain a telluric 
absorption spectrum at an airmass difference of < 0.1 
relative to the science target. The science spectra were 
corrected for telluric absorption using the IDL program 
xtellcor (Vacca et al. 2003), which removes model H i ab- 
sorption lines from the AOV spectra before application 
to the science data. All spectra were flux normalized af- 
ter modeling the continuum by the fitting of a low-order 
polynomial. 

Table 1 lists the 2MASS near-infrared positions and 

photometry of the new WR stars, and the mid-infrared 
photometry from the Spiteer/GLIMPSE database. Ta- 
ble 2 lists the instruments used to obtain the spectra and 
the dates of observation. 

4. ANALYSIS AND WR CLASSIFICATION 

Figures 3 9 show the NIR spectra of the new WRs; 
Figures 3 6 present the medium-resolution TSPEC and 
SpeX spectra, and Figures 7-9 present the lower- 
resolution OSIRIS spectra. Owing to the relatively low 
signal-to-noise ratio (S/N) of the J-band OSIRIS data, 
we present only H and K spectra for the stars observed 
with this instrument. WR subtypes were assigned follow- 
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ing the criteria of Figer ct al. (1997) and Crowthcr ct al. 
(2006), who demonstrated that WN subtypes can be clas- 
sified quantitatively via the equivalent width (EW) ra- 
tios of the following helium emission lines: He II (A1.012 
fj,m)/Re I (A 1.083 iim) in the Y band, and He ii (A2.189 
/xm)/He u+Btj (A 2.165 fim) and He ii (A2.189 iim)/Re 
I (A 2.112 fim) in the i^-band. For WN subtypes of WN6 
or earlier, the emission of high-ionization N v, relative to 
He i/C iii/N II, can also be used as a subtype diagnostic 
(Crowther et al. 2006). WC stars and their subtypes are 
classified in a manner similar to that of the WN stars via 
the EW ratios of the following carbon emission lines: C 

III (A0.971 ^m)/C II (A 0.990 /xm) in the Y band, and C 

IV (A2.076 ^m)/C ill (A 2.110 /xm) in the K band. For 
WN stars, decreasing subtype is associated with an in- 
crease in the strength of He II emission relative to He I; 
while for WC stars, decreasing subtype is associated with 
an increase in the strength of C iv emission relative to C 
III, and C III relative to C ii. We measured EW values 
using the SPLOT routine in IRAF, whereby the individ- 
ual flux contributions of blended emission lines may be 
recovered if the central wavelengths of the emission arc 
known. However, emission-line "deblending" was only 
necessary for the extraction of carbon line strengths in 
the isT-band spectra of the WC stars that exhibit very 
strong emission lines of C IV and C III. Tables 3 and 4 
list the EW values for prominent emission lines in the 
spectra of WN and WC stars, respectively, and the re- 
sulting subtype classifications. 

The new sample contains a total 38 WN stars and 22 
WC stars. The majority of WN stars have subtypes 
in the WN6 7 range, exhibiting He ii (A2.189 /xm)/Hc 
Ii-|-Br7 (A2.165 /im) ratios in the range of «0.4~1.2. Ear- 
lier WN stars, such as MDM23 and MDM52, exhibit 
slightly higher ratios of 1.6 and 2.5, respectively, consis- 
tent with earlier subtypes in the range of WN4-5. For 
the latest WN subtypes. He ii is substantially weaker 
than for the earlier WN types, as indicated by a compar- 
ison of their relatively low He II (A2.189 /im)/He Il+Br7 
(A2.165 /xm) ratios of <0.2. 

Stars MDM8 and MDM19 have very weak features of 
He II in their spectra. In the JC-band, MDM8 appears 
to have this transition in absorption, which may be re- 
garded as a feature more indicative of an O supcrgiant. 
However, in the i/-band He ii possibly appears as a weak 
P Cygni profile at A1.692 /xm, which would warrant the 
WN classification. For star MDM19, He ii at A2.189 /im 
appears to have a definite absorption component, while 
there might be some indication of an emission compo- 
nent having a P Cygni profile, although it is too near the 
noise level to be regarded as a secure detection. Thus, the 
spectral properties of MDM8 and MDM19 are mutually 
consistent with hot emission-line supergiants, possibly in 
transition to late WN stars. As an additional diagnostic, 
the strength of the Br7 emission may help discriminate 
between O and WN types. Spectroscopic analysis of O 
and hydrogen-rich WNh and WNha stars in Bohannan 
& Crowther (1999) suggest that O supergiants, such as 
HD 151804, have EW (Br7) < 10 A, whereas WN9ha 
and WN9h stars exhibit higher strengths of «20 A and 
w50 A, respectively. Following this rationale, we tenta- 
tively classify MDM8 as a WN9h star and MDM19 as a 
WN9ha star, albeit with caution. Higher S/N data with 



increased spectral resolution is necessary to assign these 
stellar classifications with more confidence. 

Star MDM50 has a particularly unique spectrum, 
which is dominated by emission lines of helium that are 
relatively narrow, compared with the other WRs. He i 
at A2.058 /xm is the strongest fine, with FWHM«300 km 
s~^, just above the 250 km s^^ velocity resolution of the 
R « 1200 spectrum. The lower FWHM values of the 
other lines in the spectrum are consistent with the veloc- 
ity resolution limit of the data, which suggests that the 
lines could be nebular in origin. In the if -band, two re- 
solved lines are apparent near A2.165 /im; their centroids 
lie at AA2.1622, 2.1654 /im, which suggests that they 
could both be He transitions (see Figer et al. 1997, their 
Table 2), specifically He i 7^L-4^D at A2.1622 /xm, and 
He I 7^D-4}F at A2.1653 /xm, although we identify these 
particular transitions with caution. There is also signif- 
icant He II emission at A2.189 /xm and A1.012 /im. The 
7?-band spectrum of MDM50 appears relatively devoid of 
any hydrogen emission features, which might also be an 
indication that the ii'-band emission feature at A2.1654 
/xm is a transition of He and not Br7. There is, how- 
ever, a potential emission line of H i (n=10-4) at A1.736 
/xm. There are also weak, but significant, features of low- 
ionization emission lines of Na i at A2.206 /xm and Mg 

II at A2.138 and A2.144 /im, implying that the emission 
at A1.736 /im could also have contribution from Fe ii at 
the same wavelength (e.g., sec Morris et al. 1996). The 
presence of low-ionization features, accompanied by high- 
ionization He II, is somewhat puzzling. The diagnostic 
line ratios we have used for the other WRs in our sample 
would imply a WN8 spectral type for MDM50. How- 
ever, the overall character of its spectrum and the nar- 
rowness of its spectral lines are not characteristic of any 
known WN8 star. The near-infrared spectrum is more 
similar to WR122 (aka NaStl), classified as Ofpe/WN9 
star by Figer et al. (1997), but re-classified as a hotter 
WR cloaked in an ionized nebula by Crowther & Smith 
(1999). The spectrum of MDM50 is inconsistent with 
either of these phenomena. Ofpe/WN9 stars typically 
exhibit much stronger emission from H I, particularly for 
the Paschen series in the H band (Morris et al. 1996). No 
H I emission has been detected from MDM50. Moreover, 
the cloaked WR NaStl exhibits a very strong thermal ex- 
cess at near- and mid-infrared wavelengths; MDM50 does 
not. Finally, as we will demonstrate in the next section, 
the photometry for this source is also inconsistent with 
the supcrgiant luminosity class. 

4.1. WC Candidates for Thermal Dust Emission 

Late-type WC stars are commonly associated with 

thermal emission from hot circumstellar dust (T ^ 500- 
1000 K), which supplements the free- free emission pro- 
duced by the stellar wind (Cohen et al. 1975; Williams 
et al. 1987). The significant near-infrared continuum 
excess can reduce the strengths of emission lines, partic- 
ularly in the K band. Stars having C ill (0.97 /xm)/ C 

III (2.11 /im) ratios that are relatively high (SS8-10) may 
be regarded as candidates for hot dust emission (WCd 
stars), since there will be more dilution in the K band 
than in the Y band. If Y band data are not available, 
stars that have generally weak lines in the K band, as 
quantified by the sum of He i (2.06 /xm) C iv(2.08/xm)-|-C 
iii(2.11/xm)-|-He I (2.06 /xm), can also be regarded as can- 
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didates for emission-line dilution by dust. However, stel- 
lar continuum emission from an unresolved OB compan- 
ion can also dilute emission lines of the WR. Thus, the 
sum of i^-band line strengths is not a definitive indi- 
cator of dust dilution; the G III (0.97 fim)/ C III (2.11 
^m) ratio is a much more reliable diagnostic. Seven of 
the WC8--9 stars in our sample have the K-band cover- 
age necessary for the calculation of the C ill (0.97 yum)/ 
C III (2.11 ^m) ratio. MDM40, MDM49, and MDM60 
have relatively high values of 9.8, 8.9, and 8.7. Of those, 
MDM40 and MDM60 also appear to be outliers having 
additional excess in the (Kg — [8.0]) versus (J — Ks) 
color space shown in Figure 1. There is also marginal 
evidence of this additional excess in the (H — Ks) ver- 
sus (J - H) color space. MDM40 and MDM60 exhibit 
relatively weak emission lines in the K band, shown in 
Figure 6, potentially the result of dilution from the ther- 
mal dust excess. MDM49 exhibits very strong lines in 
the K band, however. MDM18, MDM24, MDM26, and 
MDM58 also have relatively weak lines, yet, examining 
their colors with reference to Figure 1, these stars do not 
exhibit excesses with strength comparable to confirmed 
WCd stars; MDM40 and MDM60 do, however. In con- 
clusion, MDM40 and MDM60 are the only newly iden- 
tified WC stars to exhibit strong evidence for thermal 
dust emission, so, we classify them as WC8-9d stars. 
For HDM18, HDM24, HDM26, HDM49, and HDM58, 
the evidence for hot-dust emission is less convincing, so, 
we classify them as WCd? stars. 

5. EXTINCTION AND DISTANCE 

Using the photometric data from Table 1 and the 
spectral type classifications in Tables 3 and 4, we com- 
puted the interstellar reddening of each source by adopt- 
ing intrinsic (J — Ks)o, {H — Ks)o, and Mk values for 
WR subtypes at known distances from Crowther et al. 
(2006). We calculated the average values for (J — Ks)o, 
{H—Ks)o, and Mx for four WN7 stars that reside within 
the cluster Westcrlund 1: Stars B, D, G, and P, obtaining 
(J - Ks)o = 0.37, {H - Ks)o = 0.09, and Mk = -5.18 
mag. The same calculation was performed for the WN6 
stars Q and U in that cluster, yielding (J — Ks)o = 0.37, 
{H - Ks)o = 0.11, and Mk = -4.65 mag. We ap- 
plied these values to calciilatc the extinction and distance 
values of the WN6 and WN7 stars in our sample. For 
all other WR subtypes we adopted the calibration pho- 
tometry values from Crowther et al. (2006; their Table 
Al). Table 5 includes the adopted color and Mk^ val- 
ues, which were applied to the observed photometry to 
calculate the color excesses Ej^k^ a-^d Hh-k^ caused 
by interstellar extinction. Employing the extinction re- 
lation of Indebetouw et al. (2005), two values of A^^ 
were calculated, given by 

A-Ks = 1-821o;23Sh_K3 

and 

where E{J — Kg) = ( J — JCs)obs — [J — ^s)o) and simi- 
larly for the H — Kg colors. The two results were aver- 
aged to obtain a final Ak, estimate. Assuming values of 
Mk, based upon WR subtype;, we calculated the distance 
modulus, DM, and the physical distance, D, to each WR 
star by using the following expressions 

DM = i^, - Mk, - Ak, 



and 

D (pc) = lo(°^+5)/^ 

Assuming a distance to the Galactic center of i?o = 8 
kpc (Reid 1993), we also calculated the galactocentric 
distance of each star, given by 

Rg^ = rI + D^ cos^(6) - 2Ro D cos{l) cos{b), 

where I and b are the Galactic longitude and latitude, 
respectively. The results of these calculations are listed 
in Table 5. The uncertainties are significant. Follow- 
ing a rationale similar to that made by Sliara ct al. 
(2009), we attribute the bulk of the distance errors to 
the adopted average values of Mk^ from Crowther et 
al. (2006), which, although assigned to given subtype 
groups, exhibit significant scatter within those groups 
(«0.5 1 mag, depending on subtype). There is addi- 
tional uncertainty as a result of scatter in the adopted 
intrinsic colors (which is related to the intrinsic scatter 
in absolute magnitude for each subtype bin, and is a well 
known problem for hot stars with strong winds) and in 
the adopted extinction law. It can be argued that the ex- 
tinction properties along any highly reddened individual 
sightline are not well represented by a mean extinction 
curve; however, the differences in the redundant calcu- 
lations of Ak, , derived separately for the J — Kg and 
H — Kg colors, can be used as an indicator of the accu- 
racy of our results. The average difference between these 
two values of Ak, is «0.14 mag for WN stars and «0.23 
mag for WC stars, which is significantly less than the po- 
tential scatter in the adopted values of Mk, ■ Based on 
these considerations, the distance errors are, in general, 
estimated to be ^25%, although they might be as high as 
R:!40% in some cases. The WC9 stars, in particular, may 
have the largest spread in absolute photometry if any of 
these stars are accompanied by a hot-dust excess. How- 
ever, most of the WC8-9 stars that we have found lie in 
the WR sweet spot. Dusty WC stars should exhibit an 
additional c;olor shift in these diagrams. Therefore, it is 
reasonably safe to assume that many of the WC8-9 stars 
we have found can be regarded as not significantly dusty. 
Wc adopt Mk, = —6.3 mag for non-dusty WC9 stars, 
which is consistent with the measured values for Star E 
in the cluster Westcrlund 1 (Crowther ct al. 2006b) and 
qF 76 in the Quintuplet cluster (Figer et al. 1999). StiU, 
we present the distances for all WC9 stars with caution. 

The derived extinction and distance for MDM50 is in- 
consistent with its preliminary classification WR. 
Although the spectral morphology and line ratios suggest 
a late WN spectral type, the large distance of ^12 kpc 
that results from assigning the associated luminosity is 
very suspect, as it is grossly inconsistent with the pres- 
ence of the bright optical counterpart (B = 14.2 mag) 
and the low value of extinction (Ak, = 0.37 mag) de- 
rived for this source. Instead, the relatively blue color 
suggests that this object lies much closer to Earth and, 
as such, is intrinsically much less luminous than a late- 
type WN star. We will resume discussion of this object 
in Section 7.2. 

6. X-RAY DETECTIONS 

Nine of the WRs presented in this work are associated 
with X-ray point sources found in the Chandra archive 
(Evans et al. 2010). Table 6 lists the X-ray sources and 
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their broadband fluxes. Four sources lie within the field 
containing the young clusters Danks 1 and 2 (Danks et 
al. 1984; Georgelin et al. 1988); this field also contains 
the known colliding-wind binary WR48a (Zhekov et al. 
2011). The WN stars HDM9 from Hadfield et al. (2007) 
and WR19a (van der Hucht et al. 2001) were also iden- 
tified in the Chandra archive, and their basic X-ray data 
are presented in Table 6 as well. 

WN-type stars typically exhibit an admixture of emis- 
sion components from cool {kT < 1 keV) and hot (fcT > 
2 keV) thermal plasmas (Skinner et al. 2010). The cool 
plasma is predicted by wind shock models (Lucy & White 
1980), while the hotter component requires wind collision 
in a binary (Luo ct al. 1990) or some unidentified mech- 
anism. WC stars in binaries are typically X-ray sources, 
while in solitude they are almost never detected, prob- 
ably as a result of high opacity of WC winds to softer 
X-rays (Oskinova et al. 2003; Skinner et al. 2006, 2009). 
With the exception of the brightest colliding-wind bi- 
naries, the derived X-ray luminosities for WN stars ap- 
proximately scale with the stellar bolometric luminosity 
(Lx ^ 10~^iboi) and, perhaps more importantly, with 
wind momentum (Skinner et al. 2010); the establish- 
ment of this relationship suggests that intrinsic X-ray 
generation is a ubiquitous feature of WN stars. Bona 
fide colliding-wind binaries, however, commonly exhibit 
X-ray luminosities above the trend, sometimes by wl-2 
orders of magnitude, and typically exhibit a harder X- 
ray spectral component (e.g., see Zhekov et al. 2011, 
Anderson et al. 2011). 

To estimate the X-ray luminosities for the WR X-ray 
sources of our sample, we first converted the derived 
value of Ak^ from Table 5 to an equivalent hydrogen ab- 
sorption column, using the relation from Cardelli, Clay- 
ton, & Mathis (1989; Ay = 8.9Ak), and the relation 
from Predehl & Schmitt (1995; Nj^/Ay = 1.79 x lO^^ 
cm~^ mag~^). We then assumed that the X-ray emission 
from these WRs is thermal, having an average tempera- 
ture somewhere between 0.6 and 2 keV, consistent with 
range exhibited by WRs (e.g., see Skinner et al. 2010 and 
references therein). Using the PIMMS^ software package, 
we then calculated the unabsorbed 0.5 7 keV X-ray flux 
for the WRs and derived the resulting luminosity by ap- 
plying our derived physical distances from Table 5. The 
results are presented in Table 6. Most of the WRs have 
X-ray luminosities in the range of ^ 10^^—10^^ erg 
s^^, which is typical (Mauerhan et al. 2010b, Skinner et 
al. 2010). 

We also derived the X-ray luminosity for WR48a. ob- 
taining Lx = (3-19) X 10'^^ erg s~^. This is compara- 
ble to the value independently obtained by Zhekov et al. 
(2011), who modeled the X-ray spectrum of this object 
and derived = (8 34) x lO^** erg s^^ from a fit to a 2- 
componcnt thermal plasma having kT^^^, = 1 kcV («60% 
of the unabsorbed flux) and A:Thot = 3 keV (~40% of the 
unabsorbed flux). Their model result implies that the av- 
erage plasma temperature is indeed between 1-2 keV, as 
we have assumed. The mutual consistency implies that 
our simple method of approximating X-ray luminosities 
is sound. 

An X-ray counterpart to the WN7h star WR19a 
(van der Hucht 2001) was also identified in our cross- 

^ http://heasarc.gsfc.nasa.gov/Tools/w3pimms.html 



correlation of WR candidates with the Chandra archive. 
Using the 2MASS photometry (J = 9.08, H = 8.13, 
Kg = 7.50 mag) and intrinsic photometry from Crowther 
et al. (2006; [J - Ks]o = 0.13, [H - Ks]o = 0.11, 
Mk^ = —5.92 mag), we calculated an extinction of 
Apc, = 0.96 mag {Ay = 8.6 mag) and a spectropho- 
tometric distance of 3.1 kpc, which implies — (2- 
10) X 10'^^ erg s^^. Adopting a bolometric correction of 
BCk^ — —4.4 mag, which is the implied average value 
for two WN7-8h stars in Martins et al. (2008), we ob- 
tained Lboi = 10^'° Lq and Log(Lx/iboi) between —6.6 
and —7.3. Shara et al. (1991) measured a visual ex- 
tinction of Ay = 8.1 mag for WR19a, similar to our 
value, but calculated much larger possible distances of 
11.2 and 15.2 kpc, with their uncertainty being mainly 
the result of the My and B — V colors that they adopted. 
Applying their distances docs not have a significant ef- 
fect on the resulting Lx/iboi ratio, but implies an X-ray 
luminosity that is s::^16 25 times higher than our value 
in Table 6, which is still within the range of X-ray lu- 
minosities exhibited by WRs. Thus, the X-ray proper- 
ties do not effectively constrain the distance to WR19a. 
Still, the discrepancy between our distance to WR19a 
and that of Shara et al. (1991) is somewhat puzzling. 
To derive our distance value, we adopted Mk^ = —5.92 
mag from Crowther et al. (2006), which is the average 
value of four WN7-9h stars that exhibit individual M^^ 
values between —4.89 and —6.83 mag. If either of the 
distances derived by Shara et al. (1991) are correct, it 
would imply = —8.71 or —9.38 mag; both values 

are anomalously luminous for a WR in the infrared. So, 
cither WR19a is unusually luminous for a WN7h star, 
or the distances in Shara et al. (1991) are significantly 
overestimated. 

7. DISCUSSION 

7.1. WR Distribution and Association 

Figure 10 illustrates the galactocentric distribution of 
all known WRs with respect to the Galactic model of 
Churchwell et al. (2009); Figure 11 illustrates the lon- 
gitudinal distribution of these WRs; and. Figure 12 il- 
lustrates their galactocentric radial distribution. Table 
7 lists the current census of the Galactic WRs included 
in these figures. The Galactic distribution in Figure 10 
is likely to be affected by significant scatter, which is 
certainly a partial consequence of the distance uncer- 
tainties. Nonetheless, an overall concentration toward 
major spiral arms can be seen, while several regions of 
enhanced WR density also appear to punctuate the lon- 
gitudinal distribution. However, we stress caution while 
interpreting Figures 10-12, since the coverage of WR sur- 
veys is not evenly distributed across the Galactic plane. 
Still, it is noteworthy that, although we (i.e., Hadfield 
ct al. 2007, Mauerhan et al. 2009, and this work) have 
sampled more of our WR candidates from the side of 
the Galactic plane that is north of the Galactic center 
{I = 10°-65°), we have experienced a slightly, but sig- 
nificantly, higher WR return rate for the southern side 
of the Galactic center. W; attribute this to the Galactic 
spiral arm distribution with respect to our line of sight 
through the Galactic plane. For example, the Scutum- 
Centaurus Arm, one of the two main Galactic arms il- 
lustrated in Figure 10, and the Sagittarius Arm, both 
run approximately parallel to our line of sight between 
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/ « 300° 330° and I « 28G°-295°, respectively. As a 
result of this, there is more spiral arm coverage per unit 
longitude on the southern side of the Galactic center than 
on the northern side, at least when considering the near- 
half of the Galaxy, to which our survey has mainly been 
limited. Since massive stars are presumably concentrated 
in spiral arms, wc should expect a higher observed sur- 
face density of WRs on the southern side of the Galaxy 
than the northern side. Thus, although we approach the 
result with caution, owing to the non-uniformity and in- 
completeness of WR surveys across the Galactic plane, 
wc arc confident that the peaks of enhanced WR surface 
density in the longitudinal distribution of Figure 11 are 
manifestations of real physical structure in the Galaxy. 

The largest WR concentration is apparently toward 
the direction of the starburst cluster Westerlund 1, near 
I « 335° 340°; this line of sight pierces the Scutum- 
Centaurus Arm perpendicularly, and beyond that, runs 
approximately parallel to the Norma Arm of the Galaxy, 
which contains the f«24 WRs within Westerlund 1 and 
wlO more-distant, field WRs that lie w3-4 kpc from the 
Galactic center. We experienced our highest WR re- 
turn rate in this region of longitude; during a particu- 
larly remarkable interval of observing time, while sur- 
veying WR candidates from the WR color "sweet spot," 
near I « 340°, we confirmed 6 WRs consecutively, with- 
out interruption from non-WR "contaminants." Based 
on this experience, we regard our remaining wlO unob- 
served WR candidates from the sweet spot that lie near 
this longitude as particular promising. 

Although not selected from our color-based survey, 
there is a strong concentration of WRs toward the Galac- 
tic center. The Central Molecular Zone, which com- 
prises the central half-kiloparsec of the Galaxy, contains 
92 known WRs (see Mauerhan et al. 2010a and refer- 
ences therein), «70% of which are contained within the 
three known starburst clusters in the region: the Arches, 
Quintuplet, and Central clusters. The region as a whole 
contains ~20% of the known Galactic WR population. 
For this region of the Galaxy, Spitzer and 2MASS ob- 
servations suffer from extreme stellar confusion and very 
high extinction {Ay > 30 mag), which often prevents the 
detection of a J-band counterpart in 2MASS. The result- 
ing degradation of point-source photometry has limited 
our color-based selection to identify only one of the 92 
WRs known in this region. Thus, our method of WR 
selection is currently not effective in the Galactic center 
or fields of comparable stellar confusion and extinction, 
such as, much of GLIMPSE-II. 

For the Galactic plane on the northern side of the 
Galactic center there are also several noteworthy concen- 
trations of WRs. There is a peak near I « 30° in Figure 
11. This line of sight pierces the base of the Scutum- 
Centaurus Arm, where it intersects the Galaxy's Long 
Bar, and also contains the W43 starburst region. Vigor- 
ous star formation has occurred in this part of the Galaxy 
over the last ~20 Myr, as evidenced by the presence of 
the WR-bearing cluster W43 (Blum et al 1999), four clus- 
ters of red supergiants between I = 24° 29° (Figcr et 
al. 2006; Davies et al. 2007; Alexander et al. 2009; Clark 
et al. 2009; Negueruela et al. 2010), and an additional 
WR-bearing cluster at {l,b) = (28?46,0?32; Mauerhan 
et al. 2010c). Thus, the relatively high concentration of 
WRs here is not surprising. 



Figures 13-16 show the 2MASS /r,,-band field images 
for new WRs that do not appear to be associated with a 
local (<2') stellar density enhancement (exceptions will 
be discussed in Section 7). The relative "isolation" of the 
new WRs motivates the question as to whether all WRs 
form exclusively in clusters and associations or form via 
an alternate mode of relatively isolated birth as well. It is 
possible that some of the WRs originated within clusters 
and associations, but were subsequently ejected dynam- 
ically, either by inter-cluster gravitational encounters or 
following the supernova explosion of an erstwhile com- 
panion (Moffat & Isscrstedt 1980; Moffat ct al. 1998). 
Indeed, it has been noted that approximately half of all 
O stars exhibit kinematical evidence for runaway veloc- 
ities (de Wit et al. 2004, 2005, and references therein; 
Dray et al. 2005), and the same might be true for some of 
the WRs presented here, since they are descended from 
O stars. 

We must also consider the possibility that stellar clus- 
ters of relatively low density and mass, when compared 
to massive, starburst systems, such as Westerlund 1, can 
persist undetected as a density enhancement because of 
background confusion. It is likely that many elusive clus- 
ters could populate the Milky Way, while being unde- 
tectable within existing near-infrared surveys, such as 
2MASS. Such clusters have been referred to as "false 
negative" clusters (Popcscu & Hanson 2009; Hanson et 
al. 2010b) In the following subsection, we present the 
discovery of WRs in one such false negative cluster that 
was not found by any previous cluster search algorithms. 

7.1.1. A New Stellar Cluster at G 37. 5 1-0 46 

Figure 17 shows the if^-band field image of the WC8 
stars MDM53 and MDM54, which are separated by «7". 
Their respective distances of 6.2 and 6.0 kpc in Ta- 
ble 5 are consistent with equivalent distance, given the 
« ±0.5 mag uncertainty in the adopted value of Mk, 
that we noted in Section 5, which results in a « ±25% 
uncertainty in distance. The average distance of f«6.1 
kpc, implies a physical separation of «0.2 pc. The 
stars lie within a slight stellar density enhancement at 
(Z,6) = (37?51, — 0?46) that we noticed in field acqui- 
sition images with TSPEC while observing at Palomar 
Observatory. This possible cluster was not selected by 
any previous searches (i.e., Dutra & Bica 2001; Mercer et 
al. 2005), and so the system qualifies as a false- negative 
cluster with respect to those searches. Motivated by 
the reasonable suspicion that these WC8 stars occupy a 
massive stellar cluster, we spontaneously acquired near- 
infrared spectra of the brightest neighboring field stars, 
which are also marked in Figure 17. Table 8 lists the 
near-infrared photometry, estimated spectral types, ex- 
tinction values, and distances for the surrounding stars. 
Figure 18 presents their if-band spectra. 

Stars A, B, D, F, and G exhibit iiT-band spectral char- 
acteristics that are very similar to the B0I-B3I stars in 
cluster Sgr 1806-20 (Bibby et al. 2008), which have 
{H — Ks)o — —0.08. Adopting this value for the stars 
in cluster G37.51— 0.46, we derived extinction values for 
stars A, B, D, F, and G, and present them in Table 8. 
The average of the spectroscopic parallactic distances of 
the neighboring WC8 stars MDM53 and MDM54 (6.1 
kpc) implies Mk^ values in the luminosity class of I-II 
for stars A, B, D, F, and G. 
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Stars C, E, and H exhibit the CO bandhead character- 
istic of late-type M stars. Their derived extinction val- 
ues are significantly below those of the early-type stars, 
which implies that they lie in the foreground and are un- 
related to the system. Thus, they are not likely to be 
cool supergiants within the cluster. Even if we were to 
assuming intrinsic colors of MOIa stars, which are the 
bluest subtypes of the MI class, the extinction values 
still remain significantly lower than the values for the 
OB stars. Thus, the available data suggests that the 
late-type stars C, E, and H lie in the cluster foreground. 

The distance of 6.1 kpc places the G37.51— 0.46 cluster 
within the massive-star formation region near the inter- 
section of the Scutum-Ccntaurus Arm and Galaxy's Long 
Bar (Churchwell et al. 2009), where we noted earlier the 
enhancement in the WR surface density and the presence 
of numerous massive stellar clusters. Additional spec- 
troscopy of surrounding stars and photometric analysis 
of the field will be required to determine the full extent of 
the evolved massive-star population in this false-negative 
cluster. 

7.1.2. New WRs Near Banks 1 and 2 and Implications 

A relatively high concentration of WR candidates were 
found in the vicinity of the G305 H II region, which con- 
tains the known stellar clusters Danks 1 and 2 (Danks 
ct al. 1984; Gcorgehn ct al. 1988), shown in Fig- 
ure 19. A spectacular threc-lobcd wind-blown bubble, 
known to surround the clusters, exhibits evidence for 
vigorous, triggered massive star formation in the form 
of H2O, OH, and methanol masers that are coincident 
with numerous, embedded, infrared sources (Clark & 
Porter 2004). Several WRs arc already known to exist 
within the G305 complex, including the coUiding-wind 
binary WR48a (WC6-I-0B) and another WC8 star iden- 
tified by Mauerhan et al. (2009). We obtained spectra 
of 7 additional candidates during May 2010, shown in 
Figure 19. Stars MDM4 and MDM6 have WC8 and 
WC7 spectral types, respectively, while stars MDM3, 
MDM5, MDM7, and MDM8 have spectral types in the 
WN8-9 range. Only one of the candidates we observed 
in this region was revealed to be a Be "contaminant." 
Four additional candidates remain unobserved in the 
field covered by Figure 19, yet they are significantly 
fainter than the confirmed WRs in this region, which 
is somewhat inconsistent with them being WRs within 
the G305 complex. If they are members of this sys- 
tem, it is, perhaps, more likely that they are Be stars. 
They have the 2MASS designations J13120279-6243224 
{Ks = 12.6 mag), J13121640-6239543 {K, = 12.4 
mag), J1321223498-621456 [K, = 11.3 mag), and 
J13125377-6238308 {Ks = 11.9 mag). Of these 
unidentified WR candidates, we suspect that 2MASS 
J1321223498-621456 has the highest hkclihood of being 
a WR, since it is the brightest of the remaining candi- 
dates and closest to the main cluster center in Figure 
19. 

The Danks 1 and 2 clusters lie at an estimated distance 
of KiZ.b kpc, behind an extinction screen of Ax ~ 1 mag 
(Bica et al. 2004) . Our derived extinction and distances 
for the newly confirmed stars in Figure 19 average out 
to Ax^ = 1.2 ± 0.3 mag and D = 3.6 ± 1.1 kpc, consis- 
tent with the previously published values for the cluster. 
However, several of the individual values difii'er signifi- 



cantly from the average. Our derived distance of 6.9 kpc 
for the WC7 star MDM6 is particularly deviant, having 
nearly twice the accepted value for the clusters; yet, it 
has an extinction value that is consistent with the value 
for the region. Thus, either MDM6 lies at a distance sig- 
nificantly beyond the clusters and is not physically asso- 
ciated with them, or the adopted value of Mx = —4.59 
mag from Crowther et al. (2006) used for WC5-7 stars 
is unsuitably bright for this particular star. In fact, as- 
suming that the infrared colors that we used are appli- 
cable, MDM6 is expected to have Mx = —3.1 mag if it 
lies at the accepted cluster distance of 3.5 kpc (Clark & 
Porter 2004). The discrepancy cannot be the result of 
unaccounted-for dust emission, since the additional red- 
dening and Ks-h&nd fiux that are associated with hot 
dust emission would lead to an underestimation of dis- 
tance, not an ovcrestimation. The same would be true if 
there were a bright, unresolved companion contributing 
to the flux. Thus, we are forced to consider the possi- 
bility that WC6 stars exhibit a broader range of Mx^ 
values than we have anticipated. If we exclude MDM6 
from the calculation of average distance to G305, then 
we obtain a value of 2.9±0.8 kpc (« 28%), which im- 
plies that our anticipated, average distance uncertainty 
of «30% using the spectroscopic parallax method is a 
reasonable estimate. 

Including the WC6+0B binary WR48a, the results of 
our survey bring the total number of known WRs associ- 
ated with the G305 H 11 region to eight objects. Clark & 
Porter (2004) suggested that the Danks 1 and 2 clusters 
have succeeded in clearing away the majority of their na- 
tal material, implying that significant photon leakage is 
taking place in tandem with the ionization of the sur- 
rounding three-lobed bubble. This suggestion, coupled 
with the integrated radio fiux emanating from the region, 
led those authors to the conclusion that a large popula- 
tion of massive stars having a total Lyman continuum 
(Lyc) photon output greater than the equivalent of 31 
07V stars is likely to populate this environment, supply- 
ing at least 3.2 x 10^° Lyc s~^. The new WRs presented 
here, in addition to massive stars presented in Clark et 
al. (2011), certainly represent a significant fraction of the 
population presumed to exist by Clark & Porter (2004). 
In fact, by adopting N^y^ values of WR subtypes from 
Figer et al. (1999), we approximate that the eight WRs 
shown in Figure 19 supply 3.3 x 10^° Lyc photons s~^ 
(averaging over subtypes from WRs in Figer ct al. 1999, 
we find that log N^,, w 49.8, 49.7, and 49.5 for WC6 -7, 
WC8, and WN8 9 stars, respectively). Therefore, the 
ionizing flux of G305 inferred by Clark & Porter (2004) 
could now be fully accounted for. Of course, owing to 
the natural photon leakage out of the G305 complex and 
the obvious presence of candidate supergiant stars that 
remain unidentified within the clusters, there is undoubt- 
edly more Lyc fiux being generated in this region than 
we have currently estimated. 

It is interesting to find that five of the WRs associ- 
ated with G305 are relatively well separated from the 
centers of the Danks 1 and 2 clusters. MDM3, 4, 5, 6, 
and WR48a, as a group, appear to be centered on Danks 
1, and are separated from the cluster by «20", imply- 
ing a physical separation of f«20 pc (assuming the clus- 
ter distance of 3.5 kpc). The distribution of these stars 
is rather puzzling. As the most highly evolved massive 
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stars in the G305 system, these WRs presumably have 
the most massive progenitors, whose orbits should have 
relaxed near the cluster center as a result of dynamical 
friction. However, within the dense cluster centers there 
is a relatively high probability of gravitational encounters 
between other massive stars and binaries, which leads to 
an increased probability of orbital "heating" or cluster 
ejection. This could explain the extended distribution of 
MDM3, MDM4, MDM5, MDM6, and WR48a. We spec- 
ulated earlier that such processes provide one potential 
explanation as to why so many of the new WRs pre- 
sented in this work appear to be "isolated" objects, pos- 
sibly runaways from their birthplaces in clusters (Moffat 
& Isserstedt 1980; Moffat et al. 1998). The WRs sur- 
rounding the Danks clusters are particularly interesting 
in this regard, as they could be objects in the early phase 
of running away from their birth sites. 

7.1.3. A New WN9 Star Near Cluster [DBS2003] 179 

The WN8-9 star MDM32 is apparently associated 
with the stellar cluster [DBS2003] 179 (Dutra & Bica 
2003), which is separated from the star by «16". The 
field shown in Figure 15 includes stars MDM32 (near 
field center) and the cluster (lower- right of field center). 
Borissova et al. (2008) identified seven O stars in the 
cluster, having luminosity classes ranging from dwarfs to 
supergiants, in addition to an Ofpe/WN9 star (Object 4 
from their paper). Those authors estimated a distance 
of 7.9 kpc and an average extinction vahie of Ay = 16.6 
mag, which corresponds to Ax = 1.9 mag, assuming the 
Cardelli, Clayton, & Mathis (1989) law. Although our 
extinction value of Ak = 2.0 mag for MDM32 is consis- 
tent with the value estimated for the cluster by Borissova 
et al. (2008), our measured distance oi D = 6.4 kpc is 
significantly lower than theirs, which they measured by 
averaging spectroscopic parallax distances for the O su- 
pergiants in the cluster. If their distance of 7.9 ± 1.1 
kpc is the correct value for MDM32, then it should have 
Mk S3 —7.3 mag, which is notably higher than than the 
adopted value of Mk = —6.87 mag for the brightest late 
WN stars from Crowther et al. (2006). This discrepancy 
is not surprising, given the potential wO.S-l mag uncer- 
tainty in Mk^ that we discussed in Section 5. The dis- 
tances to the Ofpe/WN9 star Object 4 from Borissova et 
al. (2008) and MDM32 are likely to be similar; although 
saturated in H and K, the J-band and mid-infrared pho- 
tometry for Object 4 is very similar to MDM32, which 
is evidence favoring mutually close proximity to cluster 
[DBS20003] 179. In conclusion, we physically associate 
MDM32 with this cluster, and our favor leans toward 
the distance of 6.4 kpc that we estimated for the sys- 
tem. However, more observations will be necessary to 
establish a more accurate distance. 

7.1.4. Identification of a WN6 Star in Cluster 

[MFD2008] 

The WN6 star MDM34 hes within the cluster 
[MFD2008], discovered by Messineo et al. (2008). who 
associated the system with the high-energy gamma-ray 
source HESS J1813— 178. Those authors recently con- 
ducted a more thorough spectroscopic survey of candi- 
date cluster members (Messineo et al. 2011), identifying 
MDM34 (their star #7), in addition to almost two dozen 



other hot massive stars, and a red supergiant. They clas- 
sified MDM34 as a WN7 star, although the value EW 
(2.189 iim)/ EW (2.166 /zm)=1.5 and the significant de- 
tection of N V emission at A2.10 fjbva. favor the earlier 
WN6 subtype that we have assigned. Our measured ex- 
tinction of =0.88 mag and distance of 3.1±0.8 kpc 
to MDM34 is consistent with their extinction value and 
average spectro-photometric distance {Ak, = 0.83 mag; 
3.6 ± 0.7). Messineo et al. (2011) also derive an X-ray 
luminosity of 1.3 x 10^^ erg s~^ for MDM34, which is 
consistent with our estimated range of (2-22) x 10^-^ erg 
s~^. 

7.2. The uncertain nature of MDM50 

The photometry of MDM50 implies that this star's 
light suffers relatively low extinction and, thus, lies at a 
relatively small distance, suggesting that MDM50 is sig- 
nificantly less luminous than a typical WR or evolved su- 
pergiant. Furthermore, its narrow lines imply a relatively 
weak, low- velocity wind with v <250 km s~^, set by the 
resolution limit of our spectrum. Yet, the significant de- 
tection of multiple transitions of He ii presented in Figure 
3 implies high excitation, and an effective temperature 
that would place MDM50 within the O spectral class. 
Yet, the detection of weak, low-ionization features im- 
plies a separate, cooler emission component, which could 
arise from circumstellar, nebular material. MDM50's lo- 
cation in the color-color and color-magnitude plots in 
Figure 1 suggest that, although this star is one of lowest 
extinction stars in our sample, its infrared excess is con- 
sistent with the color distribution of the non-WRs (pre- 
smued Be stars) and WNs, although it lies on the blue 
edge of the broad color space. Thus, although we do not 
know the intrinsic color of MDM50, it is probably rea- 
sonable to presume that its intrinsic excess is similar to 
the bluest WRs in our sample, which have J — Kg 0.2 
mag and H — Kg ~ 0.1 mag. If so, then applying the 
method of extinction calculation from Section 6 yields 
Ak^ = 0.38 mag, using both the J — Kg and H — Kg col- 
ors. Assigning Mk values ranging from O dwarfs (« —4 
mag) to supergiants (~ —5.5 mag) we obtain possible 
distances of 3.3 kpc and 6.5 kpc, respectively. Exam- 
ining the distances of eight WRs in Table 5 that have 
Ak, < 1 mag shows that they have an average distance 
of 2.7 kpc with a standard deviation of 0.8 kpc. There- 
fore, the photometry is inconsistent with the supergiant 
luminosity class, since an extinction of Ak, ~ 0.4 mag 
would be anomalously low for a distance of 6.5 kpc. The 
Mk « — 4 mag is more appropriate, and is consistent 
with a mid-to-late OV star [L ~ 10^ Lq). 

Moving forward with the 0-star hypothesis, MDM50 
is similar in K-ham\ spectral morphology to the variable 
06V:pe star HD 39680 (note that the "pe" qualifier is 
He I emission), which is dominated by He I emission at 
A2.058 /um, while also exhibiting weak He ii at A2.189 
(Um (Hanson et al. 1996). HD 39680 is a low-amplitude 
variable star, exhibiting small fluctuations of amplitude 
0.13 mag (Lefevre et al. 2009). MDM50 is also pho- 
tometrically variable. Eight observations of calibration 
fields for the Sloan Digital Sky Survey spaced throughout 
1 year revealed MDM50 to be a variable star exhibiting 
variations of wO. 1-0.2 mag at R band (Henden & Stone 
1998), similar to the variability amplitude of HD 39680. 
It is, therefore, plausible that MDM50 is also an OV:pe 
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star like HD 39680. If so, the reason for selection of this 
star by our color criteria is unclear. It may be the result 
of excess emission from an ionized, free- free emitting cir- 
cumstellar disk similar to those surrounding the Be class 
(Negueruela et al. 2004). It is plausible that the low- 
ionization features of Na i and Mg ii that we detected 
in the MDM50's spectrum arc generated in such a disk, 
especially since we do not expect these features to arise 
from the same emission component as the high-excitation 
He II lines. As there is currently no firm evidence for 
the existence of ionized disks around Oe stars (Vink et 
al. 2009), the confirmation of such a structure around 
MDM50 would be particularly interesting. 

We also noted in Section 4 the spectral similarity be- 
tween MDM50 and NaStl, the latter of which has been 
classified as WR cloaked in a massive r] Car-like nebula 
(Crowthcr & Smith 1999). NaStl also exhibits relatively 
narrow, nebular emission lines, including He ii emission. 
However, the mid- infrared SED of NaStl exhibits a very 
large thermal excess from dust, while MDM50 exhibits a 
far weaker excess component, which is inconsistent with 
the presence of a massive dusty envelope. 

The spectrum of MDM50 shares some similarities with 
the spectra of Luminous Blue Variable (LBV) stars in 
their hot state (e.g., see Clark et al. 2011 and references 
therein). However, the absence of significant H I emission 
features in our infrared spectra is somewhat inconsistent 
with a LBV classification. Still, given that MDM50 is 
a confirmed variable emission-line star, its connection to 
the LBV phenomenon is worthy of speculation, and we 
should not rule out this possibility. If it is a quiescent 
LBV, it would have to be one of the least luminous such 
stars known, having L ^ 10^ Lq. 

In conclusion, although we regard the Oe classifica- 
tion for MDM50 as most appropriate, given the available 
data, the true nature of this star remains uncertain, and 
it deserves more detailed spectroscopic investigation and 
photometric monitoring. 

8. SUMMARY AND CONCLUSIONS 

We have spectroscopically detected 60 new Galactic 
WRs using the method of infrared color selection, and 
also by cross-correlation with X-ray point source cata- 
logs. Our WR detection rate of «20%, using the broad- 
color zones defined in Figure 1, is a significant improve- 
ment over previous efforts, which had WR detection rates 
below 10% (Hadfield et al. 2007; Mauerhan et al. 2009). 
The improvement is attributable to our refinement of the 
color selection criterion, including new color-magnitude 
restrictions, and the inclusion of X-ray emission, which is 
not a ubiquitous property of the Be "contaminants" also 
selected by our color method. Including the Be stars, it is 
remarkable that the broad color space returns early-type 
emission-line stars almost entirely, at a rate of w95%. By 
selecting WR candidates from a smaller color space of the 
"WR sweet spot," which contains significantly fewer Be 
stars, we experienced a detection rate of >50%. For the 
X-ray emitting candidates, we detected 10 WR sources 
after observing 23 objects, implying a detection rate of 
«40%. We presented their X-ray properties, which are 
consistent with known, solitary WRs, as well as those in 
some colliding- wind binaries. 

Most of the new WRs appear "isolated" in the field; 
some could be runaways, others might have formed in 



relative isolation, or could be members of "false nega- 
tive" clusters that are difficult to distinguish against the 
dense stellar background of the Galactic plane. We have 
discovered at least one such cluster, containing two WC8 
stars and several OB supergiants. We also identified sev- 
eral WRs in the known stellar clusters Danks 1 and 2, 
[DBS2003] 179, and [MFD2008]. The resuhs imply that 
our method of selecting candidate WRs from the field 
could provide the most effective means of detecting WR- 
bearing false-negative clusters in the Galaxy. 

The latest results bring the total number of known 
Galactic WRs to 476, which is «7-8% of the total ex- 
pected population of the Galaxy, estimated to be as 
high as W6000-6500 objects (Shara et al. 2009). The 
current WN/WC ratio of the Galactic sample is »1.5. 
Although WR subtype ratios are commonly referenced 
in tests of stellar evolutionary models (e.g., Meynet & 
Maeder 2005), they can only have meaning for complete 
or subtype-unbiased samples of WRs. There may be rea- 
son to suspect a subtype bias in the current sample. For 
example, even though the earliest WNE types (WN3-4) 
have dense winds (strong free- free excesses and broad He 
I and He II lines) , they are intrinsically the faintest of the 
WRs, by 1-2 mag (Crowther et al. 2006), we would ex- 
pect surveys of the inner Galaxy to have more difficulty 
recovering these objects relative to, say, WN7-9 stars. 
Thus, the Galactic longitudinal and radial distributions 
of WRs as a function of subtype, shown in Figures lib 
and 12, are simply the observed distributions, and should 
not be misconstrued or over-interpreted. On the other 
hand, the observed surface density of all known WRs 
with respect to Galactic longitude (see Figure 11a) does 
appear to reflect true physical structures of the Galaxy, 
including prominent massive-star formation sites, such 
as W43 and the Westerlund 1 region, and portions of the 
sky where major spiral arms run approximately paral- 
lel to our line of sight, e.g., the Scutum-Centaurus and 
Sagittarius arms in the range / « 280°-315°. 

There are ?a4900 unobserved WR candidates remain- 
ing in our sample from the GLIMPSE-I survey area 
(/ = 295°-350°, 10°-65°; |6| < 1°). If the overall broad- 
color detection rate of «20% persists, than we expect 
that almost 1000 WRs could potentially be detected in 
future confirmation observations. However, the observa- 
tions that we have conducted thus far have given priority 
to bright sources over fainter ones, ii's-band magnitudes 
are typically in the range of 8 11 for oTir sample, while 
average extinction values arc w2 mag. Thus, assuming a 
generic absolute Ks magnitude of w — 5 for WRs implies 
that we are mainly sensitive to WRs lying at distances 
of <8 kpc from the Sun, on our side of the Galaxy. Ex- 
panding the survey volume to include WRs all the way 
to the Solar Circle on the far side of the Galaxy would 
require sampling candidates having Ks ^ 13 mag. How- 
ever, we presume that as these fainter candidates are tar- 
geted, the associated degradation in both their 2MASS 
and GLIMPSE photometric qualities and the increase in 
stellar confusion will result in a significantly decreased 
WR confirmation rate. A decline in photometric quality 
is also expected for the inner GLIMPSE-II survey area 
(1^1 < 10° from the Galactic center), owing to the rela- 
tively high source confusion in the innermost regions of 
the Galaxy. There are ~1700 candidates in this area, 
although the WR detection rate has not yet been de- 
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termincd here, since the carhest color-based surveys by 
Hadficld ct aL (2007) failed to detect any WRs in the 
GLIMPSE-II fields. Continued surveys in this region are 
thus necessary in order to estimate the number of WRs 
that are contained within the GLIMPSE-II point-source 
catalog and the extent of the Galactic volume that our 
survey is sensitive to. 

The color-based search for Galactic WRs is currently 
the most efficient method for identifying these objects 
in the Galaxy. Although the observational prescription 
described in this work has resulted in the greatest WR 
confirmation rate produced thus far, the infrared phe- 
nomenology behind the color selection method deserves 
further study, in order to establish an improved under- 



standing of the WR subtype selection biases inherent in 
this technique. 

This publication makes use of data products from 
the Two-Micron All-Sky Survey (2MASS), which is a 
joint project of the University of Massachusetts and 
IPAC/Caltech, funded by NASA and the NSF. The re- 
search was based on observations made with the Spitzer 
Space Observatory, the Hale 5 m Telescope at Palo- 
mar Observatory, Cerro Tololo International Observa- 
tory, and the NASA Infrared Telescope Facility (IRTF). 
We thank the referee for a very insightful report, which 
helped to improve the manuscript. 
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13 
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66 
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86 
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65 


50 


18521900- 
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283 
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-0 


18649 


32 


87 


-0 


28 


9 


75 


9 


29 


8 


98 


8 


69 


8 


47 


8 


23 


7 


93 


51 


18524369+0008415 


283 


18205 





14495 


33 


21 


-0 


22 


12 


49 


10 


68 


9 


61 


8 


63 


8 


19 


7 


99 


7 


68 



USNO opt. 
(mag) 



Bl Rl 
counterpart mag 



0273-0506079 
0272-0487929 

0272- 0487952 

0273- 0506154 



0422- 



0440- 



0651194 



0524226 



20 



15 



20 



0721- 


-0782909 


18 


63 


16 


57 


0771- 


-0547089 


20 


65 


17 


43 


0769- 


-0553041 


15 


86 






0796- 


-0404402 


17 


24 


14 


41 


0798- 


-0402571 


16 


79 






0842- 


-0402447 


18 


36 


16 


96 


0849- 


-0402164 


19 


18 


15 


71 


0872- 


-0545092 


17 


96 


14 


30 


0898- 


-0397750 


14 


21 


12 


32 



13 



62 



05 



TABLE 1 — Continued 



Star 


Sonrc p 


R.A. 


Decl. 


I 


f, 


J 




s 


f3 61 


\A 51 


f5 81 


fS 01 

[O.UJ 


USNO ont 


Bl Rl 


TV /n~\A /r 
MUM 


(2MAbD J) 


(deg, J2000) 


(deg) 


(deg) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


counterpart mag 


52 


18540312+0124508 


283.51297 


1.41415 


34.49 


0.06 


12.81 


11.06 


9.99 


9.02 


8.62 


8.37 


8.09 






53 


19012661+0351553 


285.36087 


3.86544 


37.51 


— 0.46 


14.70 


12.67 


11.25 


10.05 


9.51 


9.27 


9.06 






54 


19012711+0351542 


285.36295 


3.86514 


37.51 


-0.46 


14.25 


12.27 


10.99 


9.77 


9.34 


9.09 


8.88 






55 


19083809+0928210 


287.15875 


9.47248 


43.31 


0.53 


12.22 


10.81 


9.94 


9.04 


8.64 


8.45 


8.16 






56 


19161838+1246493 


289.07663 


12.78035 


47.12 


0.39 


12.95 


11.82 


10.77 


10.22 


9.72 


9.50 


9.24 






57 


19174121+1129190 


289.42182 


11.48859 


46.13 


-0.51 


15.30 


13.00 


11.53 


10.53 


9.91 


9.59 


9.33 






58 


19202932+1412061 


290.12221 


14.20167 


48.85 


0.16 


14.34 


11.98 


10.63 


9.46 


8.93 


8.65 


8.35 






59 


19225361+1408498 


290.72343 


14.14711 


49.07 


-0.38 


12.67 


10.91 


9.68 


8.82 


8.23 


8.04 


7.73 






60 


19225446+1411280 


290.72697 


14.19109 


49.11 


-0.36 


12.83 


11.02 


9.58 


8.05 


7.61 


7.20 


7.02 






61 


19300530+1746010 


292.52211 


17.76694 


53.08 


-0.18 


12.71 


11.07 


10.12 


9.14 


8.72 


8.45 


8.13 







Note. — Photometric data is from the 2MASS (JHK^) and Spitzer/ GLIMPSE (3.6-8.0 //m) catalogs. MDM38 was independently discovered by Motch et al. (2010). MDM34 was independently 
identified in Messineo et al. (2011). 
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TABLE 2 

Spectroscopic Observations of WRs 



Star Observation Date Telescope/Instrument 
MDM (UT) 



1 


2010-05-23 


23: 


;09 


SOAR/OSIRIS 


2 


2010-05-24 


00: 


;02 


SOAR/OSIRIS 


3 


2010-05-27 


00: 


;16 


SOAR/OSIRIS 


4 


2010-05-24 


22: 


:38 


SOAR/OSIRIS 


5 


2010-05-24 


23: 


;28 


SOAR/OSIRIS 


6 


2010-05-24 


23: 


;39 


SOAR/OSIRIS 


7 


2010-05-25 


00: 


;00 


SOAR/OSIRIS 


8 


2010-05-25 


00; 


;13 


SOAR/OSIRIS 


9 


2010-05-23 


00; 


;27 


SOAR/OSIRIS 


10 


2010-05-24 


01; 


:51 


SOAR/OSIRIS 


11 


2010-05-23 


01; 


;03 


SOAR/OSIRIS 


12 


2010-05-23 


03; 


;44 


SOAR/OSIRIS 


13 


2010-05-25 


01: 


;04 


SOAR/OSIRIS 


14 


2010-05-25 


01: 


;52 


SOAR/OSIRIS 


15 


2010-05-26 


03: 


;32 


SOAR/OSIRIS 


16 


2010-05-23 


06: 


;06 


SOAR/OSIRIS 


17 


2010-05-26 


03: 


;14 


SOAR/OSIRIS 


18 


2010-05-26 


03: 


;52 


SOAR/OSIRIS 


19 


2010-05-26 


04: 


;10 


SOAR/OSIRIS 


20 


2010-05-26 


04: 


;28 


SOAR/OSIRIS 


21 


2010-05-27 


04: 


;38 


SOAR/OSIRIS 


22 


2010-05-27 


05: 


;59 


SOAR/OSIRIS 


23 


2010-05-27 


05: 


:24 


SOAR/OSIRIS 


24 


2010-05-27 


05: 


;40 


SOAR/OSIRIS 


25 


2010-05-26 


09: 


;05 


SOAR/OSIRIS 


26 


2010-05-26 


09; 


;18 


SOAR/OSIRIS 


27 


2010-05-27 


07; 


;38 


SOAR/OSIRIS 


28 


2010-05-27 


08; 


;15 


SOAR/OSIRIS 


29 


2010-05-27 


09; 


;01 


SOAR/OSIRIS 


30 


2010-05-26 


05; 


;11 


SOAR/OSIRIS 


31 


2010-05-27 


09; 


;24 


SOAR/OSIRIS 


32 


2010-05-25 


06: 


;21 


SOAR/OSIRIS 


33 


2009-08-06 


07: 


;42 


IRTF/SPEX 


34 


2009-08-05 


06: 


;15 


IRTF/SPEX 


35 


2009-07-13 


08: 


;42 


HALE/TSPEC 


36 


2009-07-13 


08: 


;25 


HALE/TSPEC 


37 


2009-07-14 


06: 


;03 


HALE/TSPEC 


38 


2009-07-14 


05: 


;13 


HALE/TSPEC 


39 


2009-07-15 


04: 


;52 


HALE/TSPEC 


40 


2009-07-15 


09: 


;38 


HALE/TSPEC 


41 


2009-08-06 


07: 


;10 


IRTF/SPEX 


42 


2009-07-15 


05: 


:14 


HALE/TSPEC 


43 


2009-07-14 


06: 


;50 


HALE/TSPEC 


44 


2009-07-15 


09: 


;16 


HALE/TSPEC 


45 


2009-08-04 


05; 


;41 


IRTF/SPEX 


46 


2009-07-14 


08; 


;52 


HALE/TSPEC 


47 


2009-08-06 


06; 


;12 


IRTF/SPEX 


48 


2009-08-06 


05; 


:40 


IRTF/SPEX 


49 


2009-08-06 


06; 


:53 


IRTF/SPEX 


50 


2009-08-05 


09; 


;22 


IRTF/SPEX 


51 


2009-07-15 


10: 


;40 


HALE/TSPEC 


52 


2009-07-15 


07: 


;10 


HALE/TSPEC 


53 


2010-06-21 


10: 


;23 


HALE/TSPEC 


54 


2010-06-21 


10: 


;42 


HALE/TSPEC 


55 


2009-07-15 


06: 


;32 


HALE/TSPEC 


56 


2009-07-15 


06: 


;51 


HALE/TSPEC 


57 


2010-06-21 


04: 


;50 


HALE/TSPEC 


58 


2010-06-21 


04: 


;37 


HALE/TSPEC 


59 


2009-08-06 


09: 


;30 


IRTF/SPEX 


60 


2009-07-15 


11: 


;20 


HALE/TSPEC 


61 


2009-07-14 


10: 


;50 


HALE/TSPEC 



Note. — Facilities have the following spectral resolutions (A/(5A): SOAR/OSIRISaiSOO; IRTF/SPEXRil200; HALE/TSPEC«2600. 



TABLE 3 

Equivalent Widths of Prominent Emission Lines for WN Stars 



05 



Star 
MDM 


He II 
1.012 


He I 
1.083 


He 11+P7 
1.093 


He II 
1.163 


He il+P^ 
1.281 


He II 
1.476 


He II 
1.692 


He I 
1.700 


He I 
2.058 


N V 
2.110 


He i/Niil 
2.115 


He Il+Br7 
2.165 


He II 
2.189 


1.012/ 
1.083 


2.189/ 

2.165 


2.189/ 
2.115 


Type 


1 










30 












3 


18 


5 




0.28 (WNS) 


1.67 (WN7) 


WN7-8 


2 










42 


6 


3 


9 


2 




11 


37 


6 




0.16 (WNS) 


0.55 (WNS) 


WNS 


3 










45 


5 


2 


9 


3 




20 


42 


4 




0.10 (WNS-9) 


0.20 (WNS-9) 


WN8-9 


5 










15 


5 


2 


17 


8 




22 


54 


4 




0.07 (WN9) 


0.18 (WNS-9) 


WN9 


7 










19 


2 


<1 


2 


<1 




9 


17 


-1 




0.05 (WN9) 


0.1 (WN8-9) 


WN9 


8 










55 


<1 


<1 


8 


6 




10 


62 


<1 




< 0.02 (WN> 9) 


< 0.1 (WN>9) 


WJN9n 


n 

y 










00 


ii 


n 

y 


ol 


1 n 
iU 




OK 
zo 


QQ 
00 


li 




0.13 (WNS) 


0.44 (WNS) 


WiNO 


1 n 














1 Q 


4U 


1 A 
14 




( 


AP. 
40 


Ofi 

ZD 




0.57 (WN6-7) 


0.71 (WN7) 


WiN / 


1 A 
±4 












1 no 


<5 ( 


9^; 
zo 


9n 

zy 





r^A 
04 


/in 
4y 


1 nf; 
iUD 




2.16 (WN3-5b) 


1.96 (WN5-7) 


\sn\TKK 
WIN Ob 














zoo 


J.UZ 


blend 


9n 

ZU 


1 A 


K7 
1 


04 


yo 




1.11 (WN6— 7b) 


1.63 (WNS— 7) 


VV iMOb 


16 










51 


31 


g 


44 


17 




48 


54 


29 




\J ,\J^ l_ W IN D — t ) 


\J-\J ^ W IN t — 0) 


WN7 


19 
















2 


1.2 




6 


26 


0.9 




0.03 (WN9) 


0.15 (WNS— 9) 


WNgha"* 


20 












27? 


28 


(blend)-)- 


4 




19 


32 


18 




0.56 (WN6— 7) 


0.95 (WN7) 


WN7 


21 










42 


51 


■(—30 


^blpnHl— i 


29 


20 


24 


37 


69 




X.OU (WJNO— 7d) 


\J.^I (WN4— D) 


WN5h'' 

VV IN OL> 


23 












115 


■(—86 


( blend J— > 


9/) 

z4 


16 


47 


26 


229 




8.81 (WN4b) 


4.87 (WN4-6) 


WiN4— OD 


OK 

zo 












iiy 


no 

■( — yz 


(^Dlenu ) — f 


1 7 
i / 


■( — 00 


\ Dlena ) — t 




1 1 Q 

lio 




1.44 ('VVN5-7b) 


2.05 (WN5-6) 


WiNO— DD 












85 


DO 


■(—88 


( blend j—>- 


C A 

54 


6 


66 


75 


RA 

64 




0.85 (WN6-7) 


0.97 (WN7) 


W1N6 


on 
zy 










44 


1 Q 


J in 

A — iy 


(^blenu ) — f 


1 n 
iU 




OQ 
ZO 


QK 
00 


Q 






U.U<5 (WN9) 


0.13 (WNS-9) 


wiNy 


30 










15 


24 


14 


7 


4 


1 


14 


17 


32 




1.88 (WN6-7) 


2.29 (WN5-6) 


WN6'' 


31 










46 


80 


32 


23 


26 


5 


28 


41 


90 




2.20 (WNS) 


3.21 (WN5-6) 


WN5'' 


32 










81 


47 


•(-40 


(blend)-^ 


182 




18 


40 


3 




0.08 (WN9) 


0.17 (WNS-9) 


WN9 


33 


8 


544 


40 


15 


54 


14 


5 


61 


62 




58 


70 


18 


0.02 (WN9) 


0.25 (WNS) 


0.31 (WNS) 


WN8-9 


34 


81 


145 


10 


61 


27 


29 


13 


21 


6 


2 


24 


29 


43 


0.56 (WN6-7) 


1.48 (WN6) 


1.79 (WN6-7) 


WN6 


35 


64 


516 


65 


40 


53 


17 


10 


52 


24 




59 


62 


27 


0.12 (WNs) 


0.44 (WNS) 


0.46 (WNS) 


WN8 


36 


98 


375 


24 


48 


60 


34 


14 


39 


14 




47 


55 


41 


0.26 (WN7) 


0.75 (WN6-7) 


0.87 (WN6-7) 


WN7 


37 


32 


33 


5 


21 


12 


11 


5 


8 


10 




6 


13 


15 


0.97 (WN6) 


1.15 (WN6-7) 


2.50 (WN5-6) 


WN6 


38 


58 


172 


19 


39 


46 


25 


^333 


(blend)^> 


19 




43 


47 


28 


0.34 (WN7) 


0.6 (WN6-7) 


0.65 (WN7-S) 


WN7 


41 


186 


517 


27 


111 


85 


44 


8 


49 


17 




67 


65 


46 


0.36 (WN7) 


0.71 (WN6-7) 


0.69 (WN7) 


WN7 


43 


93 


1026 


190 


63 


54 


40 


•(-46 (blend)-5> 


18 




61 


55 


44 


0.09 (WN7-8) 


0.80 (WN6-7) 


0.72 (WN7) 


WN7 


44 


11 


503 


22 


21 


68 




2 


32 


41 




40 


74 


11 


0.02 (WN9) 


0.15 (WNS) 


0.28 (WNS-9) 


WN8-9 


45 












11 


4 


15 


5 




20 


22 


17 




0.77 (WN6-7) 


0.85 (WN7-S) 


WN7 


46 


93 


494 


28 


57 


73 


40 


13 


57 


18 




57 


64 


39 


0.20 (WN7-8) 


0.61 (WN6-7) 


0.68 (WN7-S) 


WN7 


47 


150 


482 


6 


102 


59 


32 


14 


43 


13 




46 


59 


39 


0.31 (WN7) 


0.66 (WN6-7) 


0.85 (WN7-S) 


WN7 


48 


132 


570 


35 


60 


80 


29 


12 


47 


18 




42 


56 


32 


0.23 (WN7) 


0.57 (WN6-7) 


0.76 (WN7-S) 


WN7 


50 


9 


29 


2 


4 


9 


3 


4 


15 


73 




20 


16 


4 


0.31 (non-WR?) 


0.25 (non-WR?) 


0.20 (non-WR?) 


non-WR?"^ 


51 


12 


149 


22 


7 


49 


4 


3 


15 


7 




22 


55 


3 


0.08 (WNS) 


0.05 (WN9) 


0.14 (WNS-9) 


WN8-9 


52 


110 


118 


10 


86 


21 


35 


23 


10 


2 


6 


19 


21 


59 


0.93 (WN6) 


2.81 (WN5) 


3.11 (WN5-6) 


WN5-6 


55 


37 


175 


21 


19 


49 


11 


6 


18 


7 




21 


53 


15 


0.20 (WN7-8) 


0.28 (WNS) 


0.71 (WN7-S) 


WN8 


61 


110 


345 


24 


82 


63 


39 


13 


32 


15 




42 


61 


53 


0.32 (■WN7) 


0.87 (WN7) 


1.26 (WN7) 


WN7 



Note. — EW values are in units of Angstroms. Uncertainties are S3lO%. The last three columns list the EW ratios of diagnostic lines, with the implied subtype following in parentheses. The final subtype is given in the 
last column. The 1.012/1.083 and 2.189/2.165 ratio diagnostics are from Crowther et al. (2006, C06), and the 2.189/2.115 ratio diagnostic is from Figer et al. (1997, F97). Final classifications respect both criterion. Stars are 
assigned the broad designation ("b") if the FWHM of He 11 at 2.189 /limis above 130 A, following C06. 

^ Spectra consistent with both O supcrgiants and hydrogen-rich WN atars (WNh or WNha) were further classified by the overall strength of Br^- emission at A2-165 ,^m (sec text). 

^ For these particular stars, the ratio of N V (2.10 /^m)/He l/NlII (2.115 p.m) is also used as a diagnostic, following C06. MDM21, for example, has a possible subtype range of WN4-7, given the C06 and F97 criterion. However, 
the spectrum exhibits 2.100/2.115=0.8, which implies WNS, according to C06. So, WNS was assigned. 

The nature of MDM50 is uncertain; the subtype diagnostics implied by the line ratios are not considered applicable. 



TABLE 4 

Equivalent Widths op Prominent Emission Lines for WC Stars 



Star 
MDM 


C III 
0.971 


C II 
0.990 


He II 
1.012 


He I 
1.083 


C IV C III 
1.191 1.198 


He l/C II 
1.700 


C IV 
1.74 


He I 
2.058 


C IV 
2.076 


C III 
2.110 


He ll/C III 
2.165 


He II 
2.189 


0.971/ 
0.990 


2.076/ 
2.110 


0.971/ 
2.110 




Type 


4 












weak 


89 




340 


118 


10 


5 




2.9 (WC8) 




458 


WC8 


6 












40 


501 




1862 


385 




80 




4.8 (wc:7) 




2247 


WC7 


11 












8 


97 




486 


106 


20 


34 




4.6 (WC7) 




592 


WC7 


12 












22 


108 


85 


256 


149 


18 


24 




1.7 (WC8) 




490 


WC8 


13 












56 


70 


109 


161 


172 


51 


51 




0.9 (WC9) 




442 


WC9 


17 












87 


184 




556 


256 


114 


57 




2.2 (WC8) 




812 


WC8 


18 












64 


32 


54 


45 


80 


73 


19 




0.6 (WC9) 




179 (d?) 


WC9d? 


22 












77 


36 


362 


79 


135 


88 


21 




0.6 (WC9) 




576 


WC9 


24 












25 


47 


25 


82 


84 


24 


28 




1.0 (WC8-9) 




191 (d?) 


WC8-9d? 


26 












29 


11 


76 


31 


65 


34 


26 




0.5 (WC9) 




172 (d?) 


WC9d? 


28 












82 


134 




352 


195 


68 


59 




1.8 (WC8) 




486 


WC8 


39 


714 


71 


80 


537 


•(-202 (blend)-)- 


84 


109 




303 


173 


67 


60 


10 (WGS) 


1.8 (WC8) 


4.1 


476 


WC8 


40 


616 


64 


26 


282 


30 114 


25 


19 


68 


44 


63 


20 


19 


9.6 (WC9) 


0.7 (WC9) 


9.8 (d) 


175 (d?) 


WC9d 


42 


743 


58 


93 


600 


2o2 69 


100 


169 




433 


260 


65 


57 


12.8 (WC8) 


1.7 (WC8) 




693 


WC8 


49 


2460 




257 


315 


•(-390 (blend)-)- 


12 


290 




1073 


276 


95 


43 




3.9 (WC8) 


8.9 (d?) 


1349 


WC8d? 


53 










91 


212 




560 


227 


52 


57 




2.5 (WC8) 


787 


WC8 


54 












49 


116 




374 


152 


35 


42 




2.4 (WC8) 




526 


WC8 


56 


1001 


89 


138 


215 


^329 (blend)-5> 


63 


311 




1202 


422 




132 




2.8 (WC8) 


2.4 


1624 


WC8 


57 










132 78 


65 


266 




1036 


226 


105 


83 




4.6 (WC7) 




1367 


WC7 


58 








358 


26 25 


53 


38 




71 


69 


60 


29 




1.0 (WC8-9) 




200 (d?) 


WC8-9d? 


59 


897 


26 


15 


550 


155 16 


95 


203 




766 


159 


54 


81 


34.5 (WC5-7) 


4.8 (WC7) 


5.6 


925 


WC7 


60 


980 


49 


73 


324 


149 12 


37 


108 




299 


113 


21 


22 




2.6 (wcs) 


8.7 (d?) 


412 (d?) 


WC8d 



Note. — Equivalent widths are given in units of Angstroms. The uncertainties are P^10% of the given value. Columns 15 and 16 list the equivalent width ratios of selected diagnostic lines from Crowther et al. (2006), and the 
implied WC subtype follows in parentheses. The final assigned subtype is given in the last column. The ratio of C iii lines given in column 17 is used as a diagnostic to select candidates for thermal dust emission, which can 
dilute line strengths in the K-band relative to the V band. Column 18 lists SK, which is defined as the sum of the equivalent widths of He i (2.06 /jm), C iii (2.11 /L*m), and C iv (2.08 Aim) in the K band. For WCS— 9 stars, this 
value may be relatively low for thermal dust emitters, but this does not necessarily imply dust emission (see text). 



18 



TABLE 5 

Extinction and Absolute Photometry for New WRs 



Star 
MDM 


Type 




{H - Ks)o 




(niSiSi) 


(uiSiSi) 




R 
(kpc) 


Rg 


1 


WN7-8 


0.13 


0.11 


1.43 


1.36 


1.39 


-5.92 


3.7 


8.0 


2 


WN8 


0.13 


0.11 


2.12 


1.97 


2.04 


-5.92 


5.5 


7.2 


3 


WN8-9 


0.13 


0.11 


1.68 


1.60 


1.64 


-5.92 


2.4 


6.9 


4 


WC8 


0.43 


0.38 


1.04 


0.76 


0.90 


-4.65 


3.2 


6.7 


5 


WN9 


0.13 


0.11 


1.36 


1.31 


1.33 


-5.92 


2.8 


6.8 


6 


WC7 


0.62 


0.58 


1.23 


0.96 


1.10 


-4.59 


6.9 


6.9 


7 


WN9 


0.13 


0.11 


1.03 


1.04 


1.04 


-5.92 


4.0 


6.6 


8 


WN9h 


0.13 


0.11 


1.02 


0.99 


1.00 


-5.92 


2.0 


7.0 


9 


WN8 


0.13 


0.11 


2.02 


2.00 


2.01 


-5.92 


8.3 


7.4 


10 


WN7 


0.37 


0.09 


2.13 


2.21 


2.17 


-5.18 


7.9 


6.8 


11 


WC7 


0.62 


0.58 


1.83 


1.27 


1.55 


-4.59 


8.3 


2.6 


12 


WC8 


0.13 


0.11 


2.98 


2.91 


2.95 


-5.92 


6.2 


6.1 


13 


WC9 


0.23 


0.26 


2.01 


1.80 


1.91 


-6.30 


9.0 


6.8 


14 


WN5b 


0.37 


0.27 


2.55 


2.31 


2.43 


-4.77 


4.6 


4.8 


15 


WN6b 


0.37 


0.27 


2.82 


2.58 


2.70 


-4.77 


4.9 


4.6 


16 


WN7 


0.37 


0.09 


2.92 


2.89 


2.90 


-5.18 


4.3 


4.8 


17 


WC8 


0.43 


0.38 


2.56 


2.15 


2.36 


—4.65 


5.3 


4.4 


18 


WC9d? 


0.23 


0.26 


2.59 


2.30 


2.45 


-6.30 


11. 9*' 


6.2'' 


19 


WN9ha 


0.13 


0.11 


3.21 


2.89 


3.05 


-5.92 


4.9 


4.2 


20 


WN7 


0.37 


0.09 


2.57 


2.52 


2.55 


-5.18 


6.5 


3.5 


21 


WN5b 


0.37 


0.27 


1.69 


1.39 


1.54 


-4.77 


10.0 


4.3 


22 


WC8-9d? 


0.23 


0.26 


2.08 


1.95 


2.01 


-6.30 


15.3 


8.5 


23 


WN4-5b 


0.37 


0.27 


2.50 


2.30 


2.40 


-4.77 


4.5 


4.3 


24 


WC8-9d? 


0.23 


0.26 


2.62 


2.39 


2.51 


-6.30 


ll.O'' 


4.8'' 


25 


WN5-6b 


0.37 


0.27 


1.95 


1.77 


1.86 


-4.77 


4.7 


4.1 


26 


WC9d? 


0.23 


0.26 


1.95 


1.63 


1.79 


-6.30 


9.8*' 


4.0'' 


27 


WN6 


0.37 


0.11 


2.12 


2.27 


2.19 


-4.70 


6.3 


3.3 


28 


WC8 


0.43 


0.38 


1.86 


1.64 


1.75 


-4.65 


4.5 


4.2 


29 


WN9 


0.13 


0.11 


2.87 


2.66 


2.77 


-5.92 


5.2 


3.7 


30 


WN6 


0.37 


0.11 


1.10 


1.14 


1.12 


-4.70 


3.4 


4.9 


31 


WN5 


0.18 


0.16 


2.26 


2.01 


2.13 


-4.41 


3.7 


4.7 


32 


WN9 


0.18 


0.12 


2.01 


2.03 


2.02 


-6.87 


6.4 


2.2 


33 


WN8-9 


0.13 


0.11 


1.85 


1.71 


1.78 


-5.92 


10.3 


3.0 


34 


WN6-7 


0.37 


0.11 


0.85 


0.91 


0.88 


-4.70 


3.1 


5.0 


35 


WN8 


0.13 


0.11 


1.42 


1.34 


1.38 


-5.92 


5.8 


3.1 


36 


WN7 


0.37 


0.09 


0.62 


0.79 


0.70 


-5.18 


3.6 


4.7 


37 


WN6 


0.37 


0.11 


0.59 


0.67 


0.63 


-4.70 


5.0 


3.8 


38 


WN7 


0.37 


0.09 


0.73 


1.03 


0.88 


-5.18 


2.4 


5.8 


39 


WC8 


0.43 


0.38 


1.60 


1.43 


1.51 


-4.65 


3.2 


5.2 


40 


WC9d 


0.23 


0.26 


2.14 


2.13 


2.14 


-6.30 


4.9'' 


4.2'' 


41 


WN7 


0.37 


0.09 


1.80 


1.91 


1.85 


-5.18 


6.1 


3.6 


42 


WC8 


0.43 


0.38 


1.23 


1.00 


1.12 


-4.65 


4.1 


4.7 


43 


WN7 


0.37 


0.09 


0.90 


1.10 


1.00 


-5.18 


3.5 


5.2 


44 


WN8-9 


0.13 


0.11 


2.06 


1.92 


1.99 


-5.92 


5.2 


4.2 


45 


WN7 


0.37 


0.09 


2.44 


2.41 


2.42 


-5.18 


6.0 


3.9 


46 


WN7 


0.37 


0.09 


0.74 


0.97 


0.86 


-5.18 


3.7 


5.1 


47 


WN7 


0.37 


0.09 


1.82 


1.83 


1.83 


-5.18 


7.8 


4.0 


48 


WN7 


0.37 


0.09 


1.82 


1.91 


1.87 


-5.18 


6.5 


3.9 


49 


WC8d? 


0.43 


0.38 


2.04 


1.98 


2.01 


-4.65 


4.6'' 


4.7'' 


50=* 


if OVe 


0.2 


0.1 


0.38 


0.38 


0.38 


-4 


3.3 


5.5 


51 


WN8-9 


0.13 


0.11 


1.84 


1.75 


1.79 


-5.92 


5.6 


4.5 


52 


WN5-6 


0.18 


0.16 


1.77 


1.65 


1.71 


-4.41 


3.5 


5.5 


53 


WC8 


0.43 


0.38 


2.02 


1.89 


1.96 


-4.65 


6.2 


4.9 


54 


WC8 


0.43 


0.38 


1.90 


1.63 


1.76 


-4.65 


6.0 


4.9 


55 


WN8 


0.13 


0.11 


1.44 


1.39 


1.42 


-5.92 


7.7 


5.8 


56 


WC8 


0.43 


0.38 


1.18 


1.22 


1.20 


-4.65 


7.0 


6.1 


57 


WC7 


0.62 


0.58 


2.10 


1.61 


1.86 


-4.59 


7.1 


6.0 


58 


WC8-9d? 


0.23 


0.26 


2.33 


1.99 


2.16 


-6.30 


8.9 


7.1 


59 


WC7 


0.62 


0.58 


1.59 


1.18 


1.39 


-4.59 


3.8 


5.5 


60 


WC8d 


0.43 


0.38 


1.89 


1.93 


1.91 


-4.65 


2.9'' 


6.5'' 


61 


WN7 


0.37 


0.11 


1.49 


1.53 


1.51 


-4.70 


4.6 


6.4 



Note. — The values of (J — Ks)q, (H — Ks)q, Mk^ were adopted from Crowther et al. (2006b) according to WR subtype (see text). Two 
independent values of Ak^ were calculated and averaged for a final extinction estimate. The distance uncertainties are estimated at fti25%, but 
may range as high as Ki40% (see text). 

^ Although the nature of MDM50 is uncertain, the available data are most consistent with an OVe classification and the intrinsic color values 
in the table (see text). 

^ The distances to these stars should be taken with additional caution, since they are candidates for thermal dust emission. As a result, the 
adopted colors and Kg-haxid absolute photometry used to make the distance calculation might not be reliable. 
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TABLE 6 
X-RAY Detections of WRs 



Star 


X-ray Source 


Fx (0.5-7.0 kcV) 




(0.5-7.0 kcV) 


Lx (0.5-7.0 keV) 


MDM 


[CXO Jl 


[ergs s ^ cm 




[10^2 cm-2] 


[ergs s ^ cm ^ 


'] 


[ergs s" 


1 


102328.8-574629 


6.08|:I« X 10- 


-14 


2.2 


(1.8-9.3) X 10- 


13 


(3-16) X 10^2 


3 


131209.0-624326 


6.68|:i X 10- 


-14 


2.6 


(2.2-12.6) X 10- 


13 


(1-18) X 10^2 


4 


131221.2-624012 




-15 


1.4 


(2.2-8.6) X 10" 


■14 


(3-11) X 10^1 


5 


131225.4-624441 


2.39f X 10- 


-15 


2.1 


(6.9-34.8) X 10- 


IB 


(6-33) X 1030 


8 


131228.5-624143 


4.io|:ii X 10- 


-14 


1.6 


(1.0-4.4) X 10" 


13 


(5-22)xl03i 


WR48a 


131239.5-624255 


5.00t:0| X 10- 


-12 


3.3"* 


(1.8-12.8) X IQ- 


11 


(3-19)xl034 


30 


164746.0-455904 


5.48t:g^ X 10- 


-14 


1.8 


(1.5-6.7) X 10- 


13 


(2-9) X 10^2 


34 


181322.4-175350 


2.04f;|| X 10- 


-14 


1.4 


(2.0-19.0) X 10" 


■14 


(2-22) X 10^1 


36 


182606.1-130410 


l.OOi;!! X 10- 


-14 


1.1 


(2.4-8.0) X 10" 


14 


(4-12)xl03i 


38 


183116.4-100924 


6.366;6| X 10- 


-13 


1.4 


(1.5-5.9) X 10- 


12 


(1-4) X 1032 


45 


184332.5-040419 


3.70t;2ll X 10- 


-14 


3.9 


(1.5-11.8) X 10" 


13 


(7-51)xl032 


HDM9 


181922.1-160309 


1.26;;f|] X 10- 


13 


1.2 


(2.8-10.0) X 10- 


13 


(4-17)xl033 


WR19a 


101853.3-580753 


8.58^J| X 10- 


-14 


1.5 


(2.1-8.6) X 10" 


■13 


(1-10)X1032 



Note. — Basic X-ray photometry data were obtained from the Chandra point source catalog (Evans ct a. 2010). Fx is the total flux within 
the 0.5—7.0 keV energy band, and includes 90% confidence values. Nh^^ is the interstellar absorption column derived by inserting the infrared 
extinction values in Table 5 into the relations of Cardelli et al. (1989) and Predehl & Schmitt (1995). F^^'^' and Lx are the absorption-corrected 
fiuxes and luminosities, respectively, derived using the PIMMS software package. The spread in the individual F^'^ and Lx values were derived 
assuming that the X-ray sources can be approximated by a single-temperature thermal plasma having kT in the range of 0.6—2 keV, which is 
generally a safe approximation for both single WRs and those in colliding- wind binaries. We estimate that the derived Lx values are accurate to 

1 dex. The table includes X-ray data for new WRs, as well as for WR48a, HDM9 from Hadfield et al. (2007) and WR19a, the latter two of which 
were not previously recognized as X-ray sources. Stars MDMl, MDM3, MDM4, MDM5, MDM8, and WR48a lie in the field of the Banks 1 and 

2 double cluster shown in Figure 19. WR48a, a known coUiding-wind X-ray binary (WC6-|-0), is included in the table for comparison; it is 2—3 
orders of magnitude brighter in X-rays than any of the neighboring sources. 

^ Using the method in Section 6, we derived the extinction of 71/^^—2.05 mag for WR48a (WC6+0) using its observed 2MASS photometry 
(J— 8.74 , H—6.80. and 5.09 mag ) and assuming the intrinsic colors for WC5— 7 stars from Crowther et al. (2006); the same calculations were 

performed for WR19a (,7—9.08 , H=S.13, and Kg— 7.50 mag ), for which we obtained Aifg=0.96 mag and D — 3.1 kpc. The corresponding value 
of Nh^^ was calculated in the same manner as the other sources. We use the accepted average distance to Banks 1 and 2 of 3.5 kpc for WR48a. 



TABLE 7 
Total Galactic WR Census 



Reference N (WR) 



Van der Hucht (2001, 2006) 


171 WN, 10 WN/WC, 113 WC, 4 WO 


Kurtev et al. (2007) 


4 WN 


Hadfield et al. (2007) 


11 WN, 4 WC 


Hyodo et al. (2008) 


1 WC 


Messineo et al. (2009) 


2 WN, 2 WC 


Liermann et al. (2009) 


2 WC 


Mikles et al. (2006) 


1 WN 


Hanson et al. (2010a) 


1 WN 


Gvaramadzc ct al. (2009) 


1 WN 


Gvaramadzc ct al. (2010) 


1 WN 


Mauerhan et al. (2007) 


2 WN 


Mauerhan et al. (2009) 


9 WN, 3 WC 


Mauerhan et al. (2010a) 


5 WN, 6 WC 


Mauerhan ct al. (2010b) 


7 WN, 2 WC 


Mauerhan ct al. (2010c) 


2 WN, 1 WC 


Motch et al. (2010) 


1 WN 


Roman-Lopes et al. (2011) 


1 WN, 1 WC 


Shara et al. (2009) 


15(13)^ WN, 26 WC 


Wachter et al. (2010) 


5 WN, 1 WC 


Anderson et al. (2011) 


3 WN 


This work 


38 WN 22 WC 


Total 


(278 WN)+(10 WN/WC)+(184 WC)+(4 WO)=. 



Note. — The table includes the WRs cataloged in van der Hucht (2001, 2006) and other references that contain all subsequent discoveries. 

Note that some of the listed works contain references therein that credit the original discoveries made by authors unlisted in this table. 



^ 2 WN stars were co-discovered by Mauerhan et al. (2009) and Shara et al. (2009); so we only count 39 WRs from the latter author. 



TABLE 8 

Photometry of Cluster Neighbors to WC8 Stars MDM52 and MDM53 



Star 


Source 
(2MASS J) 


R.A. Decl. 
(deg, J2000) 


J 


H 






(if - K)o 






11 J. dig 


111 dig 


Spectral 
TvDe 


A 


19012504+0351392 


285.354335 


3.860890 


13.491 


11.367 


10.294 




-0.08 




2.10 


2.10 


-5.73 


BO-3 (I-II) 


B 


19012634+0351471 


285.359753 


3.863086 


12.829 


10.733 


9.630 




—0.08 




2.15 


2.15 


—6.45 


BO-3 (I-II) 


C 


19012726+0351461 


285.363620 


3.862817 


13.143 


11.074 


10.102 


0.79 


0.18 


1.51 


1.44 


1.48 




M 


D 


19012736+0351587 


285.364040 


3.866331 


13.590 


11.409 


10.252 




-0.08 




2.25 


2.25 


-5.93 


BO-3 (I-II) 


E 


19012738+0352203 


285.364104 


3.872317 


13.412 


10.871 


9.673 


0.79 


0.18 


1.98 


1.85 


1.92 




M 


F 


19012804+0351581 


285.366835 


3.866142 


12.504 


10.281 


9.115 




-0.08 




2.27 


2.27 


-7.08 


BO-3 (I-II) 


G 


19012889+0351541 


285.370405 


3.865032 


14.136 


11.977 


10.927 




-0.08 




2.06 


2.06 


-5.06 


BO-3 (I-II) 


H 


19013009+0351595 


285.375405 


3.866535 


13.101 


10.806 


9.831 


0.79 


0.18 


1.66 


1.45 


1.56 




M 



Note. — Photometric data arc from the 2MASS catalog. Extinction values were derived after adopting the intrinsic colors of BO— 3 I stars from Bibby et al. (2008). Mk^ values were then 
calculated assuming a distance 6.1 kpc, which is the average distance of the two WC8 stars MDM53 and MDM54 from Table 5. The extinction values of the late-type stars were derived after 
assuming the intrinsic colors of MOIa stars from Elias et al. (1985); these colors were adopted to demonstrate that the photometry of the late-type stars is inconsistent with them being red 
supergiants within the cluster. The data suggest, rather, that they are foreground stars. 
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Fig. 1. — Infrared color-color and color-magnitude diagrams illustrating the color selection criteria of our survey. The solid black lines define the broad color 
space, discussed in the text. WRs that lie within the GLIMPSE survey area are represented by green and red square symbols. Dusty WC stars are represented 
by red inverted triangles. Our latest discoveries (this paper) are marked with an additional, open circle. Non-WRs detected in our survey, which are mainly 
Be stars, are represented by blue filled circles. The small black points are stars randomly drawn from the GLIMPSE point-source catalog. The grey shaded 
region in the upper two plots is the WR "sweet spot," discussed in the text; the confirmation rate of bona fide WRs is relatively high within this region, at 
>50%. 
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0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

Fig. 2. — Color-color diagram similar to Figure Ic, illustrating the effect that line emission has on the infrared excess 
of WRs. The colored lines connect the observed colors (filled symbols) of selected WN and WC stars to the adjusted 
colors that result from the removal of emission lines (open symbols). WN stars MDM33 (WN9, left) and MDM41 
(WN6, right) show httle color shift as a result of emission hues, while WC stars MDM42 (WC8, left), MDM56 (WC8, 
center) and MDM59 (WC6, right), show significant color shift that mainly results from strong emission lines of C ill 
and C IV in the K band. The remaining excess relative to the photospheric emission from normal dwarfs and giants 
(small black points), after the removal of emission lines, is free-free emission from the ionized winds and, possibly, hot 
dust for some WC stars. 
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Fig. 3.— Near-infrared spectra of WN stars at i? w 1200-2500 (see text regarding MDM50). 
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Fig. 4.— Near-infrared spectra of WN stars at i? 1200-2500. 
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Fig. 5.— Near-infrared spectra of WN stars at i? 1200-2500. 
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Fig. 6.— Near-infrared spectra of WC stars at i? « 1200-2500. 



28 




Fig. 7. — H and K spectra of WN stars at i? w 800. Bold face is used for the i^-band spectrum of MDM32, so as 
not to confuse the relatively strong He emission line near A2.06 /xm from this star with the other spectra in the figure. 




H and K spectra of WN stars at i? w 800. 




Fig. 9. — H and K spectra of WC stars ak Rk, 800. Bold face is used for the -fC-band spectrum of MDM22, so as 
not to confuse the relatively strong emission line near A2.06 yum from this star with the other spectra in the figure. 
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Fig. 10. — Galactic distribution of WRs and selected WR-bearing clusters plotted over the Galactic model of 
Churchwell et al. (2009). New identifications are marked as bold symbols. The figure axes are in units of kpc, with the 
location of the Sun at the origin. Dotted lines mark several values of constant Galactic longitude (l). The GLIMPSE 
survey area encompasses | ^ |< 65° on both sides of the Galactic center. The following clusters are listed clockwise from 
left with respect to the Sun: Gas 0B5 (Negueruela 2003) ;Thc Quartet cluster (Messineo et al. 2009); G37.51-0.46 
(this work; circled); G28.46+0.32 (Mauerhan et al. 2G10c); W43 (Blum et al. 1999); GLIMPSE20 (Messineo et al. 
2009); [MFD2008] (Messineo et al. 2008); Sgr 1806-20 (Bibby et al. 2008); Quintuplet, Arches and Gentral parsec 
of the Galactic center; [DBS2003] 179 (Borissova et al. 2008); Westerlund 1 (Crowther et al. 2006); Banks 1 and 2 
(Bica et al. 2004); GLIMPSE30 (Kurtev et al. 2007); Trumpler 16 (Walborn 1973); NGG 3603 (Moffat 1983); and 
Westerlund 2 (Rauw et al. 2007). The Galactic center contains sa20% of the known Galactic WRs. 



32 



GC 



W43 region, 
SC Arm 
Meets Bar 



Orion Spur 

+ Cyg 0B2 - 



WDl region 
' -rNorma Arm 



SC Arm along 
line of sight 



1/ 



V 



Sgr Arm along 
line of sight 



150 



100 50 350 300 

Galactic Longitude (degrees) 



250 



200 




150 100 50 350 300 

Galactic Longitude (degrees) 



250 



200 



Fig. 11. — Galactic longitudinal distribution for known WRs included in Table 7. The histogram in the upper panel 
is for all WRs, including those outside of the GLIMPSE survey area, and has bins of 5° longitude. The histogram in 
the lower panel shows the subtype distributions for the same WRs, and has bins of 10° longitude. 
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Fig. 12. — Galactocentric radial distribution of all known WRs included in 7. 
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Fig. 13. — 2MASS Kg-hand field images of new WRs. Tlie numbers at tlie upper left of the individual images 
represent the names introduced in Table 1 (e.g., frame 1 is for star MDMl). 
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Fig. 14.— 2MASS K^-hand field images of new WRs (con't). 
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Fig. 15.— 2MASS JsTg-band field images of new WRs (con't). 
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Fig. 16.— 2MASS isTg-band field images of new WRs con't. Note the faint nebulosity surrounding the WC8 star 
MDM59, which is presumably wind-swept interstellar medium, as indicated by Spitzer /IRAC mid-infrared images of 
this field (not presented here). 
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RIGHT ASCENSION 

Fig. 17.— 2MASS JC^-band field image of the cluster containing the WC8 stars MDM53 and MDM54, and the 
spectroscopically-confirmed, neighboring OB supergiants listed in Table 8 (Stars A, B, D, F, and G; see spectra in 
Figure 18). Stars C, E, and H are late-type stars; their relatively low extinction values suggest that they lie in the 
foreground. 
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Fig. 18. — i^-band spectra of selected stars in the vicinity of the clustered WC8 stars MDM53 and MDM54. Stars 
A, B, D, F, and G have spectral morphologies similar to known BO-3 I~II stars from Bibby ct al. (2008). The spectra 
of stars C, E, and H exhibit features of late-type M stars; their relatively low extinction values imply that they lie in 
front of the cluster. 
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Fig. 19. — 2MASS Kg-hand image of the field containing the clusters Danks f (right) and 2 (left). The newly identified 
WRs are marked, along with the two previously identified WRs in the region. Unidentified WR candidates in the field 
are also marked, although they are significantly fainter than confirmed WRs associated with the clusters. 
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Fig. A.l. Example and K spectra of non-WR "contaminants" that have selected by our method, which are 
dominated by a series of H i emission hnes. We conclude these are Be stars, based on the similarities with the Be 
spectral atlases of Clark & Steele (2000) and Steele & Clark (2001). Be stars mimic WRs in infrared color space, 
owing to similarities between the physical parameters of the free- free emitting Be disks and WR winds (see text). 



